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ABSTRACT 


The  notatkxial  stabiltty  of  adual-qnn.  qiusi^gid,  axisymmeiric  qncecnft  oontaimng 
a  driven  rotor  is  analyzed.  The  pmpose  is  to  examine  a  revised  enosy-sinkstabili^dieofy 
that  properly  accounts  for  the  energy  contiibotioo  oi  the  motor.  An  inctxisistency  in  the 
development  ditqproves  die  existing  energy-sink  theory's  assumpticm  diat  die  motor  ci  the 
system  contributes  exacdy  enough  energy  to  t^fset  die  frictional  losses  between  die  rotor 
and  die  platform.  Using  die  concept  of  core  energy,  die  revised  stabili^  criteria  for  a  dual¬ 
spin,  quasi-rigid,  axisymmetric  qracecraft  ctxitaiiiing  a  driven  romr  is  derived.  An 
expression  for  nutation  angle  as  a  ftincdon  of  core  energy  over  time  is  then  determined. 
Numerical  simulations  are  used  to  verify  the  revised  energy-sink  stability  theory.  The  dual¬ 
spin,  quasi-rigid,  axisymmetric  system  presented  by  D.  L.  Mingori  was  chosen  for  the 
simulation.  Equations  for  angular  momentum  and  total  energy  were  necessary  to  validate 
the  numerical  simulation  and  confitm  aspects  of  die  revised  energy-sink  stability  dieory. 
These  equations  are  derived  from  die  first  principles  of  dynamics  and  are  included  in  the 
uialysis.  An  explicit  relationship  for  core  energy  as  a  function  of  time  does  not  exist 
Various  models  postulating  core  energy  are  presented  and  analyzed.  The  numerical 
simulations  of  die  conqiuted  nutatitm  angles  as  a  function  of  the  postulated  ccne  energy 
compare  weU  widi  die  actual  nutation  angles  of  the  system  to  ctmfirm  the  revised  energy- 
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L  INTRODUCTION 


A  dvoaok^ical  review  of  eariy  quioecnift  9pes,  and  the  stabUiQr  criteria  developed 
fothetn.  provides  a  suffidembadcgioiind  for  die  fimdameotalcoooqKSf^  this  thesis.  The 
revised  enei|y-sink  stability  criteiioo  is  dien  presented,  and  an  equation  for  nutation  angle 
as  a  function  of  energy  diss^Mtion  over  dine  is  developed. 

A.  SINGLE  SPIN  SATELLITES 

1.  Equations  of  Motion 

The  earliest  sateUiies  took  advantage  oi  the  hex  that  stability  could  be  achieved  via 
the  'gyroscopic  stiffness*  of  a  spinning  body.  A  pceliminaiy  dynamic  model  for  the 
satellite  could  be  achieved  with  Euler’s  equations  of  motion  for  a  rigid  body.  Eu' 
moment  equation  can  be  written  as 

M  +  (1> 

dt  dt 

In  conyonent  fenn  diis  becomes 

M\  =  ^i  +  <U2%“®3^2 

Ml  s*  ^  +  o>3  hi  —  ®i  hs  (2) 

Ml  *  ^  +  ®ill2“®2hi 

Simplificatioa  of  the  model  is  accomplished  by  assuming  that  it  is  axisymmetric,  that  the 
body-fixed  axes  coincide  with  die  principal  axes  (correct  for  simple  spacecraft  with 
1 1 «  /2).  and  that  die  body  is  in  a  torque-free  environment  (a  valid  approxinoatimi  used 
dnoughottt  this  thesis).  The  qiacecraft  is  dien  rejsesented  by  Hgure  (1)  and  the  equations 
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0  - 
0  - 

0»h»i 


(3) 


Axtsynometric  Rigid  Body 
Rgnre(l) 


The  angolir  velocity  vector,  angular  momentum  vecmr,  and  die  kinetic  energy  may  be 
eiqutsaed  leqiectivety  as 


01  *mibi’i'm2h2'«’m3b3 

(4) 

b  *  /i  mi  bi  'f/2  02  b2  013  bs 

(5) 

£  •  l(fi  af  h  »3*) 

(6) 

A  simplification  of  die  notation  can  be  made.  Let /}  <B/2«/f  and/3«/s.  Bomdietliiid 
line  of  Equation  (3)  it  can  be  seen  diat  the  angular  velocity  about  the  spin  axis  b3  is 
constant,  therefore  m3  «  a>«  .  The  transverse  angular  velod^  conqionents  interchange 
between  die  hi  and  the  62  axes  but  the  magnitude  is  cmistant,  so  one  can  let 
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»  #1  ^  ♦  Us  ki.  Ttanfion, 


«  5^-»-ni«h3 

0) 

h  ■  /|S|4>/,a>,h3 

(8) 

*  /  *  «  *  +  /,* 

(9) 

E  - 

(10) 

I  is  toiqae-fiee.  Mis  constant  and  hit  fixed  in  qMce.  Because  diere  are 
no  energy  sources  or  sinks,  E  is  also  constant  Finally,  the  above  conditums  result  in 


a>  Ibemgaconstttit 

An  ex|»ession  for  the  nutation  angle  may  now  be  developed.  The  orientation  of 

die  body  axes  widi  re::  :t  to  an  inertial  reference  frame  is  desired  to  provide  a  measure  of 
die  body's  dynamic  behavior.  Because  the  angular  momentum  vector  is  fixed  in  ^lace,  the 
nutation  angle  is  defined  as  the  angle  between  die  body-fixed  axis  about  which  spin  is 
desired  and  die  angular  momentum  vector,  and  can  be  expressed  in  one  of  the  following 
fonns 
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e 

e 


(11) 

(12) 

(13) 


fa  these  etpations  die  first  and  second  tennsconeqiond  to  die  general  case  with  spin  about 
die  bs  axis,  and  the  third  term  uses  die  sinqilified  notation  to  describe  an  axisymmetric 
body,  fa  the  fecial  case  of  qnn  about  only  one  {vincipal  axis  of  an  axisymmetric  body, 
the  angular  momentum  vector  and  die  angular  velocity  vector  will  lie  on  die  qpin  axis.  Widi 
qnn  components  on  two  or  more  princ^  axes,  the  three  vectors  will  not  be  coincident. 


aUiou^  diey  still  win  Ue  on  die  same  fdane. 
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2.  Single  Spin  Satellite  Stability 

A  tocque-firee,  uisymmetric  rigid  body  with  the  body  axes  cdnciding  with  the 
princ4>le  of  inertia  will  be  stable  about  the  axis  of  eidier  die  maxiiiuini  moment  of 

inertia  or  the  minimum  moment  of  inertia.  To  prove  dns,  one  begins  with  an  artntiaiy  rigid 
body.  The  body  is  given  the  initial  oonditioa  of  steady  angular  velocity,  tuo.  about  a 
princ^  axis,  and  is  dien  perturbed  sli^tly.  b  is  assumed  diat  the  angular  velocities  about 
the  other  axes  are  small,  and  are  ^iproximately  the  same  order  of  magnitude  as  the 
perturbation  (oh  -  a>2  *  The  system  will  be  considered  stable  if  die  perturbation  does 

not  increase  over  time.  Given  an  initial  angular  velocity  with  a  perturbation, 
N  B 

S  >  oil  bi  o>2  b2 +(too'*‘£)b3,  and  given  arbitrary  inertias  l\,  I2,  h,  Euler's 
equations  cf  motkm  can  be  written  as 

0  *  /i  +(73-/2)  <»2(«>d  +  c) 

0  *=  /2  a>2  +  (/i -h)  (o>o  +  e) 

0  *  h^*{h-h)o>\  o>2 

The  equations  are  linearized  neglecting  terms  of  magnitude  £^.  Rewriting  the  equations 
by  eliminating  the  terms  e  Q>i<  £  >nd  oii  <02  results  in 


Q>2  OId 

7l 

(15) 

a>2  *  OTO 

h 

(16) 

=  fflo  +  e  «  c  «  0 

(17) 

From  Equation  (17),  (Hie  can  conclude  that  es  constant  By  differentiating  EquaticHi  (15), 
and  using  Equation  (16)  to  eliminate  012,  one  gets 

»  0  (18) 
Similarly,  from  differentiating  Equation  (16),  and  using  Equatum  (IS)  to  eliminate  (Ui,  one 
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(19) 


If>  is  imaginary.  Oh  ando^  will  increase  without  bound  over  time,  and  the  moti(»  will  be 
unstable.  Stalnlity  is  achieved  if  y  is  real.  Tbe  first  case  occurs  when  the  maximum 
moment  of  inertia  is  about  die  sfia  (ba)  axis;  then  /a  >  /i ,  /a  >  /a.  and  (/a  -  h)  ih  -  /i )  >  0. 
In  this  case  die  inertia  ratio /«//|,  defined  as  the  inertia  about  the  ^pin  axis  over  the  inertia 
about  a  transverse  axis,  is  greater  dian  one.  The  second  case  occurs  when  the  minimum 
moment  ofinertia  is  about  the  spin  (ba)  axis;  then /a  </i  ,  h<l2,  and  (/a-/2)(/a-/i)>0 
as  well  Here  the  inertia  ratio  is  less  than  one.  In  both  of  the  above  cases  7  is  real  and  the 
motion  is  stable.  However,  if  /a  is  the  intermediate  moment  of  inertia  about  the  ^in  (ba) 
axis,dien(/a-/2)(/a'-/i)<0,  7  is  imaginaiy,  and  the  motion  is  unstable. 

The  previous  model  cannot  be  ai^lied  to  a  satellite  since  die  assun^tion  of  a  rigid 
body  cannot  be  extended  to  the  spacecraft  Structural  elasticity,  liquid  jm^Uant  slosh, 
etc.,  cause  energy  dissipation  in  an  actual  spacecraft  This  spacecraft  can  be  generalized  by 
a  quasi-rigid  body  with  an  unspecified  energy  danqier  mechanism.  A  priori,  one  can 
conclude  that  energy  in  the  above  system  will  dissqNUe  until  tbe  minimum  energy  state  is 
readied.  The  kinetic  energy  of  die  ^stemwidi^dn  about  die  prindpal  axis  widimaxinuim 
moment  oi  inertia  and  qnn  about  die  prindpal  axis  with  minimum  mcmient  of  inertia  can  be 
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witeen  ia|M9Ctivdy  as 


Imm^tboatbs 

/«fa«boutb3 


03) 


The  utgfilar  momratum  is  constant  in  die  torque-fiee  case.  Thus  die  minimum  kinetic 
eneigy  state  occurs  sidien  die  rotation  is  about  die  axis  of  die  maximum  momwit  of  inertia. 
Therefore,  a  qoasi'iigid  body  is  staUe  only  when  it  is  qnnning  about  its  mgor  axis,  with  a 
corresponding  inertia  latiodiat  is  greaw  than  one. 

A  relationship  can  be  established  betwemi  the  time  rate  of  change  of  the  nutation 
angle  and  the  energy  dissipadon  (rf  a  quaa-rigid,  axisymmetiic  body.  One  must  assume 
that  the  angular  momentum  and  moments  of  inertia  of  the  quasi-rigid  body  do  not  change 
qipreciably  from  a  conqiarable  rigid  body.  For  die  genenlized  model,  with  arbitrary 
sttd  ti.  Equations  (5)  and  (6)  are  substituted  into  Equation  (12)  to  obtain 


The  time  rate  of  change  cf  die  nutation  angle  is  determined  by  taking  the  derivadve  of  the 
above  equation 


unal 


(25) 


sin(2d)(/3-/i)h* 

udiere  die  only  rale  of  change  energy  £  is  anributed  to  the  danqnng  mechanian  and  is 
written  as  £|>iomi  to  emphasize  diis  point  Because  £i>  Miri  is  negative,  the  nutation  angle 
will  decrease  only  if  /a  is  greater  dian  /i.  This  reaffirms  die  previous  conclusitxi  that  a 
quasi-rigid  body  is  Sjnn  stabilized  only  about  die  axis  <rf  maximum  moment  of  inertia.  The 
forgoing  devdopment  is  refoired  to  as  die  energy-sink  mediod. 
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B.  DUAL  SPIN  SATELLITES 

The  logical  progressioii  from  the  single  qnn  satellite  was  to  incorporate  a  de>spun 
platfonn.  This  pennitted  the  xeplacement  of  die  low  gain  omnidirectional  antenna  with 
directional  antennas  for  communication  satellites,  and  a  more  capable  spacecraft  for 
scientific  observation.  A  simple  control  system  about  the  bs  axis  would  maintain  the 
platfonn  rotating  at  a  constant  relative  rate  with  respect  to  the  rotor  ( and  would  usually 
have  the  platform  rotate  at  the  earth's  rotational  rate).  Initially,  the  platforms  were 
sufficiently  small,  and  die  overall  dimensions  of  the  satellite  were  such  that  the  inertia  ratio 
woulo  be  greater  than  cme.  For  this  type  of  satellite,  the  previously  developed  theory 
proved  adequate.  However,  as  satellites  continued  to  grow  in  size,  the  launch  vehicle 
shroud  diameter  became  a  constraint  In  order  to  provide  the  size  spacecraft  needed  to 
satisfy  mission  requirements  and  still  fit  within  die  shroud,  a  qiacecraft  with  an  inertia  ratio 
of  less  than  one  {Is  teuu  I  It  unai  <  1)  would  need  to  be  built.  From  the  previously 
developed  theory,  a  spacecraft  with  an  inertia  ratio  of  less  than  one  was  believed  to  be 
inherently  unstable.  It  was  not  until  the  development  of  the  energy-sink  theory  for  dual¬ 
spin,  quasi-rigid,  axisymmetric  spacecraft  containing  a  driven  rotor,  developed 
simultaneously  by  V.  D.  Landtm  (unpublished  wt^)  and  A.  J.  lorillo  [Ref.  1],  that  a 
spacecraft  with  an  inertia  ratio  less  than  erne  was  ctmsidered  feasible.  Several  rigorous 
stability  analyses  using  the  equations  of  motion  for  specific  dual  spinners  have  been 
performed  by  P.  W.  Likins  [Ref.  2],  D.  L.  Mingori  [Ref.  3],  and  others,  to  validate  the 
energy-sink  theory.  The  difficulty  of  a  rigmous  analysis  is  in  accurately  modeling  all  the 
forms  of  energy  dissipation.  A  more  general  and  practical  approach  was  required  to 
determine  stability,  and  the  energy-sink  theory  proved  suitable.  The  development  of  this 
theory  is  as  follows. 
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The  qieoecnft,  shown  in  Hgure  (2),  is  isramed  to  fulfill  die  ficrilowing  conditions: 

•  Both  die  losor  and  die  idatfonn  are  ixi^inmeiric 

•  Bodi  die  rotor  ud  die  platfdnn  are  qoasi-iigid 

•  The  danqdn^  mechanisms  do  not  significandy  alter  die  energy  value,  although 
die  mechanisms  will  contribute  to  TO  energy  me 

•  No  external  tonpies  are  applied 

•  The  oidy  relative  morion  is  qiin  about  the  bs  axis 

•  The  motor,  driven  die  contrcdsystet^ii^ts  just  enough  energy  to  exacdy 

offset  die  shaft  frictional  losses,  maintaining  a  constant  relative  angular 
velocity  between  the  rotor  and  the  platform 


Dual-Spin  Quasi-Rigid  Axi^onmetric  Spacecraft 
Figure  (2) 


The  Biigiiitude  of  angular  niomailiim  and  the  kinetic  eneffy  die  dual-qnn  system 
be  expressed  re^tcctively  by  the  following  equations 

A*  *  /i*ioM/ Ok* +  +//«>/)*  ■  ♦  il.  o+Li^cd/Y  (26) 

//«)/*  ««/«)/ (27) 

where  <T  is  die  rdadve  lotatioo  rate  oi  die  rotor  with  reflect  to  die  platform.  Although  the 
equatkxis  actually  represent  a  rigid  body  syswn,  diey  are  also  qipUcable  to  die  quasi-rigid 
system  because  of  the  above  assumptions.  If  the  danqnng  mechanisms  do  make  significant 
contributions  to  the  energy  of  the  system,  thoi  Equations  (26)  and  (27)  do  not  hold,  and 
the  energy-sink  criterion  will  not  qiply.  Additionally,  the  potential  energy  of  the  system 
(for  exanqile  the  energy  stored  in  die  spring  of  a  mass-spring-dashpot  damper)  will  not 
make  a  significant  contribution  to  the  total  energy  of  die  system.  Therefore,  for  the 
remainder  of  the  diesis,  it  is  assumed  that  the  kinetic  energy  of  the  system  is  effectively  the 
total  eneigy  of  die  system  (e  s  EmJI’  Because  the  system  eiqieriences  no  external  torques, 
angular  momentum  is  conserved.  Because  the  motor  contributes  no  energy  to  the  system, 
the  time  rate  of  diange  of  die  kinetic  energy  £  is  ie|»esented  by  only  the  quasi-rigid  eneigy 
dissipation  &d  total*  and  is  negative.  Differentiating  the  above  two  equations  widi  reflect  to 
time,  <me  obtains 

0  =  l^tetai  ®I  +  {it  +  //  ®1  +  //  ®/)  (28) 

£  -  ^Dtotal  °  /|io(a/®r<^  +  /i®iaV'*'//®/®/  (29) 

Eliminating  die  common  tenn  ®i  <P|  by  combining  the  above  equations 

01,  +  Is'  tUsjjl,  w,  //  ®j 

It  total 

One  m^r  now  define  the  inertial  nutation  frequency,  Ao,  the  rotcv  nutation  frequency. 


') 


+ 1,  ®,  +  //  ©/  (30) 


D  total  ’ 


Ai> 


9 


Xt  and  die  piatfonn  mnatioo  fieqoency,  X\  as  follows 


Xo 


If  (Df  +  If*  to. 


It  utai 

X  «  Ao-ttf  > 
X*  *  Xo—tOf*  * 


(31) 

• 

(32) 

•t  total 

^  IsOh’^W-IttoJioh* 

It  total 

(33) 

Hie  nutation  angle  for  a  dual-q»n  system  is  defined  as 


0  *  cos-* 

By  inclosing  die  oonditkm 


(1^)  . 


Xo>0  (35) 

the  analysis  of  nutadonal  modmi  is  restricted  to  the  following  regimi  without  any  loss  of 
generality 

(36) 

Incoqxnating  the  nutation  frequency  tenns,  the  equation  for  eneigy  dissqiation  is  written  as 

^  total  *  ^0*^0  *  ^IfXODs  —  lf  X*  its  (37) 

Recalling  the  assunqition  diat  the  rotor  and  piatfonn  are  uncoupled  about  the  bs  axis,  we 
may  incoqxmte  Equation  (37)  into  die  reaction  torques  which  taxi  to  change  the  angular 
rates 


(38) 
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CSooMiiiBi  Bqnttioii  (3(^  (37X  «ad  (39),  one  lohws  at  the  tnn^^ 


h  foul  ®f  ^ 


(39) 


DiSta^ntiating  Equation  ( 12)  and  substituting  Equation  (39)  into  it,  die  time  rate  of  diange 
of  nutation  angle  as  a  function  of  enefgy  dissqation  tales  is 


*^i.{2^)A>U  v' 


(40) 


The  eneigy-sink  equation  for  a  single  qnn-stabilized  body  is  obtained  by  letting  A  and  X' 
become  Xq  and  +  io  become  mo/  for  *  single  body.  The  definition  of  stability 


requires  die  nutation  angle  0  to  remain  constant  or  decrease  as  a  function  of  time,  so  that  0 
is  zero  or  negative.  Because  the  factors  outside  the  parenthesis  (xi  the  right  hand  side  of 
Equatitm  (40)  are  positive,  the  stabili^  criterion  for  a  dual-qiin,  quasi-rigid,  axisymmetric 
system  becomes 


^  0 

X  X' 


One  die  ftdlowing  cases  will  guarantee  sodality 


(41) 


1) 

2) 

3) 


A>0  and  A'>0 


A  >  0 ,  A'  <  0 ,  and 

A<0,  A'>0,  and 


A  specific  example  would  be  the  model  of  a  typical  communications  satellite,  a  prolate 
dual-spinner  possessing  an  inertia  ratio  of  less  than  one.  In  general,  die  rotor  nutation 
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frequency,  expitssed  as 


h  t0Ml 


(42) 


X 


would  be  negative  because  (/«-/(i0ia{)  is  negative  and  if  0/  is  rotating  at  the 

eaxtfa's  rotation  rate.  Thus,  energy  danying  in  the  rotor,  would  be  destabilizing  while 
energy  damping  in  the  platform,  would  be  stabilizing.  It  is  frtun  this  result  that 
satellites  will  have  a  damping  mechanism  placed  on  the  platfcnn  to  improve  nutational 
stability.  Such  a  danqrer  is  called  a  nutation  damper. 
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n.  PROBLEM  ranNinoN 


A.  OVERVIEW 

The  exittmg  eaetgy'tdnk  dwoiy  relies  on  levenil  assumptions,  periuqps  die  most 
important  relating  to  the  driven  rotor.  As  previously  stated,  it  has  been  assumed  diat  die 
motor,  driven  by  the  contnd  system,  iiqnns  just  enough  energy  to  exacdy  offset  the  shaft 
fiictiooal  losses,  maintaining  a  constant  idadve  angular  velocity  between  die  rotor  and  the 
pladbnn.  In  actual  systems,  contraiy  to  dus  assuiqition,  die  motor  may  add  or  remove 
energy  ftom  the  system,  dqiending  on  die  dynamics  of  the  spacecraft.  Consequendy,  a 
revised  energy  sink  stability  dieoiy,  properiy  accounting  for  die  eneigy  contribution  of  the 
motor,  is  derived.  The  revised  theory,  based  on  the  concept  core  energy,  will  remain 
consistent  with  the  existing  energy-sink  stalxlity  criterion.  Continuing,  an  equation  for 
nmation  angle  over  time,  as  a  function  of  core  eneigy,  is  developed.  Given  a  postulated 
eneigy  dissqntion  function  modeling  die  notation  dampers,  structural  elasticity,  fuel  slosh, 
etc.,  one  can  accuiitely  predict  the  nutation  angle  bduvior.  Numerical  simulation  t^D.  L. 
Mingori's  dual-spin,  quasi-rigid,  axisymmetric  system  containing  a  driven  rotor  [Ref.  3]  is 
used  to  validate  die  revised  energy-sink  stability  theory.  The  predicted  nutation  angle, 
based  on  dus  revised  energy-sink  theory,  and  the  postulated  oiergy  dissipation  function,  is 
compared  to  die  exact  nutation  angle  (rf  the  Miitgori  system.  By  using  a  suitable  postulated 
energy  dissqwtion  function,  one  can  achieve  excellent  agreement  between  the  predicted  and 
die  exact  nutation  angle.  I.  Michael  Ross  [Ref.  4]  perftnmed  diis  analysis  on  a  dual  %)in 
qrstem  with  a  damper  on  die  platform  only.  The  remainder  of  this  thesis  will  use  die  same 
analysis,  but  on  die  h^gori  system  widi  dampers  on  both  the  rotor  and  the  platfonxL 
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B.  INCONSISTENCY  OF  THE  ENERGY  SINK  THEORY 

There  exists  a  contradiction  between  die  existing  energy-sink  theoiy  and  the  nutation 
angle  derived  from  it  This  will  provide  die  motivation  for  developing  a  revised  energy- 
sink  theory  and  an  alternative  equation  for  nutation  angle  over  time.  Rom  before,  the 
existing  energy-sink  stability  criterion  can  be  expressed  as 


(41) 


aixl  die  nutation  angle  for  the  dual-spin  system  was  defined  as 


0 


(34) 


As  previously  stated,  for  a  prolate  dual-spinner  (/« iota/ //ii0ia/<  1).  energy  dissipation  in 
the  platform  is  stabilizing  and  energy  dissipaticm  in  the  rotor  is  destabilizing.  The  angular 
velocity  of  the  platform  about  the  spin  axis,  oi/,  can  be  expressed  in  terms  of  the  angular 
momentum  and  the  kinetic  energy  of  the  system.  Combiiiing  Equation  (26)  and  Equation 
(27),  one  arrives  at 


fs  tctal  *  \Is  total  lolalift  total  ~  total) 


Substituting  this  expression  into  Equation  (34)  results  in 

0  =  cos~*(±Vc) 

where  Q  is  represented  as 


h^U-ljJataL 

*t  total 


InauL 


h  total h  total 


(43) 


(44) 


Initial  conditions  at  r «  0  will  determine  the  correct  sign,  with  continuity  considerations 
maintaining  the  sign  for  all  of  r  >  0.  Additionally,  no  external  tmques  are  applied  to  this 
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^sinii,nniliiilg  in  M  being  ooBSM  DilftmnittknofEqiiaikn  (43)  results  in 


B  m  - 


i_4.(±V3)  -  ±-L-Jut#-^ 


Vl-(±V3)*  *  sin«V5  A2|i_(Lifiis() 

'  h  uialf 


(45) 


I¥om  (he  definition  of  nutation  angle,  Equatkm  (34),  and  the  conditioa  imposed  on  it. 
Equation  (36),  the  positive  sign  must  be  chosen  in  Equatioos  (43)  and  (45).  An  impoftant 
obseivation  is  made  at  diis  time.  Choosing  the  positive  sign  will  result  in  a  positive  rate  of 
change  in  the  nutation  angle,  indicating  an  unstable  condition.  The  relative  rotor  qnn  rate  is 
an  independent  variable,  and  is  arbitrarily  selected  here  as  a  constant  value  over  time. 
Therefore,  energy  dissipation  in  a  prolate  dual-^»n  ^Mcecraft  will  make  the  nutation  angle 
increase,  regardless  of  whether  the  dissipation  is  in  the  platform  or  in  the  rotor.  This  is  not 
consistent  with  the  stability  criteritm  of  Equation  (41).  Thus,  the  existing  energy-sink 
stability  criterion  contradicts  itself. 


C.  CORE  ENERGY  AND  ENERGY  DISSIPATION 

The  existing  energy-sink  stability  criterion  does  not  properly  account  for  any  energy 
that  may  be  provided  by,  or  absorbed  by,  the  motor.  To  accurately  represent  the  system, 
die  total  rate  of  change  of  energy  mist  be  written  as 

^•^010101+^  (46) 

where  is  the  rate  of  work  due  to  the  moror,  and  may  be  either  positive  or  negative  and 
the  rate  of  change  of  energy  due  to  dissqiative  elements  can  occur  in  either  die  platform  or 
the  rotor.  Recall  that  the  kinetic  energy  of  the  dual-^in  system  was  expressed  in  Equation 
(27).  If  die  work  due  to  the  motor  torque  as  a  function  of  time  is  written  in  analytical  form, 
then  the  time  rate  of  change  of  the  energy  of  the  system  due  to  all  dissipative  elements, 

£d  uitaif  can  then  be  expressed  solely  in  terms  of  die  quasi-rigid  parameters  of  the  dual-qiin 
systenL  With  this  expressicm,  the  cmidition  that  uiai  ^  0  will  result  in  the  required 
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stability  criterion.  The  difficulty  arises  in  that  to  get  die  woik  due  to  the  motor  torque  W 
(or  W),  one  needs  to  know  the  exact  dynamks  of  the  dissipative  mechanism.  In  the 
developoKnt  die  modified  eneigy-sink  stability.  Equation  (46)  will  be  used  to  detennine 
die  expression  for  teub  *tid  ultimately  derive  the  revised  stability  criterion  and  nutation 
angle  equation. 

Additionally,  the  existing  energy-sink  theory  can  be  shown  to  be  incorrect  for  both 
the  case  of  total  energy  decreasing  and  for  the  case  of  total  energy  increasing.  For 
exanqile,  if  die  rotor  is  rigid  and  total  energy  decreases.  Equation  (41)  predicts  that  the 
system  will  be  stable,  but  Equation  (45)  predicts  diat  it  will  be  unstable.  Allowing  the  total 
energy  to  increase  would  reverse  the  conditions,  but  still  show  a  ccmtradiction  between 
Equation  (41)  and  Equation  (45). 

A  modification  of  the  energy-sink  stability  theory  and  the  associated  expression  for 
nutation  angle  is  now  presented.  The  development  of  the  theory  is  fiom  I.  Michael  Ross' 
unpublished  notes.  The  basis  the  new  theory  is  centered  on  the  core  energy  of  the 
system.  As  defined  by  Hubert  [Ref.  5] 

Gore  energy  is  the  total  energy  of  the  ^cecraft  minus  the  portion  of  the  rotor 
energy  that  is  due  to  the  relative  rotatitm  between  the  rotor  and  the  platform.  It 
is  assumed  that  the  mass,  inerti^  and  motion  of  the  damping  device  are 
sufficiendy  small  that  its  energy  is  negligible  relative  to  the  ^acecraft  core 
energy.  ‘The  damper  will  be  treated  as  an  undefined  'energy  sink'  for  the 
purposes  of  the  energy  sink  analysis. 


From  the  above  statement  one  can  define  the  core  body  as  that  body  whose  inertial 
dyruunics  are  selected  for  arutlysis.  Hubert  defines  die  platform  as  the  core  body.  The  core 
energy  is  sinqily  the  rotational  kinetic  energy  of  a  fictitious  rigid  body  that  possesses  the 
inertia  prc^ierties  of  the  entire  dual-qrinner  but  nooves  in  inertial  space  exactly  like  the  core 
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body.  l^«dMl<^)toq»oecnift  the  iriaifionn  cons  eocfgy  is  defined  as 

«  ^/|  igf^  OD}^ 'f  ^/2  ipia{  4- ifigl  013'  *  ^/rtoiaf^t^ ^Artoia/ (^*^) 

Tbe  oenier  expresskm  isfior  an  aifaitnry  dual-qnn  spococah  with  tfw  qrin  axis  about  the  bs 
axis,  and  the  right  expression  is  for  a  dual-spin,  axisymmetric  spacecraft  with  the 
simplified  notation.  Extending  this  concept  10  die  rotor,  the  rotor  core  energy  is  defined  as 

£(C  *  ^/i  loMi  w*  ■*“  2  ^2  total  ^2  *2^^  "  2  ****^  ^  2 

Necessary  to  die  develt^ment  of  die  modified  dieory  is  what  will  be  termed  the 
Sqiaratitm  Axiom.  This  is  when  analysis  is  first  petfonned  with  the  rotor  considered  rigid 
and  the  platform  quasi-rigid.  Euler's  equations  are  written  for  the  rotor,  and  through 
manipulation,  an  equation  is  derived  relating  die  torque  on  the  quasi-rigid  platform  solely  in 
terms  of  platform  variables.  Then  the  platform  is  conridered  rigid  and  the  rotor  quasi-rigid. 
An  equation  is  derived  relating  the  torque  on  the  quasi-rigid  rotor  solely  in  terms  of  rotor 
variables.  These  two  separate  equations  are  then  combined  and  applied  to  a  system  in 
which  the  rotor  and  the  platfcmn  may  both  be  quasi-rigid. 

The  case  of  a  rigid  rotor  with  a  quasi-rigki  platform  is  first  analyzed.  Because  the 
rotor  is  rigid,  the  torque  applied  by  the  motOT  to  the  rotor  is  the  net  torque  on  the  rotor  and 
is  determined  fiom  Euler's  moment  equations.  For  die  case  of  the  axisymmetric  rotor, 

Tjt  •  Ta/M  »  /*a»r  (49) 

One  can  observe  that  Tp/m  »  -  Tk/m  ,  but  Tp  ^  Tp  since  the  damping  mechanism 
contributes  additional  torques  to  tbe  quasi-rigid  platfoniL  The  rate  of  wtnk  needed  by  the 
motor  torque  to  maintain  constant  relative  motion  between  the  rotor  and  the  platfoim  is 

W  Tpa  »  IsCOft  =  If  0(0/  +  ^  (50) 

By  substituting  the  platform  core  energy.  Equation  (47),  into  die  kinetic  energy  expression. 
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Eqpittian  (27),  kinetic  eneigy  may  be  expressed  as 

EmEc'^^I,<^  +  lgah'o  (51) 

The  above  equation  is  difierentiaied  to  anive  at  the  time  late  cf  change  of  die  Idnetic  eneiigy 
of  diedual-q[nn  system 

i  +  »  tc*  +  Is  <f(^' ^  *  Is  Oh'  o  (52) 

Substituting  into  diis  equation  the  rate  of  woik  done  by  die  motor  torque.  Equation  (50), 
ooeg^ 

&  ^  ic'^^  +  IsOh'6  (53) 

Conqiaring  this  to  Equation  (46),  it  can  be  seen  that 


^Dtaud^  +IsOh'  O  (54) 

Taking  the  derivative  of  Equation  (47)  to  get  the  time  rate  of  change  of  the  platform  core 
energy 


^C^ItioudOH^'^IsuiatOh'^'  (55) 

and  then  substituting  this  into  Equation  (54),  one  arrives  at  die  total  energy  dissipation  of  a 
rigid  rotor,  quasi-rigid  platform  system 


total  —  It  total  OH^'^Ig  total  Oh'  +  Is  Oh'  O 

~  It  total  Oh  Is  Oh'  +  Is'  Oh'  (Os'  +  Is  Oh' 6 


(56) 


Because  the  system  has  no  external  forces  applied,  it  remains  toique-firee.  Thus,  Equation 
(28)  can  be  used  to  eliminate  die  term  uid  arrive  at 


'D  total 


{is^  +  Ii 


total  ”  Is  total)  (^'  —  Is^ 


It  total 


(57) 
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Netfa^tlMitfi8ptalfoWiiHMtkmlflwqiiency,Bqtttik«(33),<^ 

r  -  Ao -  OV'  - * {f$tcmi-Itioimt)oh'-*'h  o 
'  A  total  /  A  total 

dial  Equation  (S7)  can  be  written  as 

total  “  (^*  + V 

Refening  to  Equation  (34),  die  nutadon  angle  can  be  written  as 


(58) 


(59) 


cos  6 


(60) 


Taking  the  derivative  and  conqiaring  it  to  die  tate  (tf  total  oiergy  diss^ation,  Equation  (59), 
die  fdlowing  relationship  can  be  established 


V  N  /  A'W 


(61) 


Because  the  rotor  is  rigid,  all  energy  dissipation  will  occur  in  die  platform,  such  diat 

—  »  |hj  Bsin  0 
X' 


(62) 


Equations  (61)  and  (62)  can  be  rewritten  by  including  the  motor  torque.  Equation  (49),  as 

(63) 


6n 

-WCsind  =  r*  +  //tti!,' -  ® 


A' 


Because  Tpim  «  -  Tg/M  (action  -  reaction  pair),  die  final  relationship  is  written  as 

Tp/M  *  (64) 

A' 

This  equation  describes  the  motor  tmque  <m  the  quasi-rigid  platform  as  a  function  of 
platfiamvaiiaUesonly.  It  can  be  seen  at  diis  point  that  the  classical  analysis  can  be 
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acilieved  by  assimiiig  no  nqiie  is  qiplied  by  the  motor,  itsiibing  m 

(65) 

r 

This  analysis  can  now  be  perfonned  on  die  qrstem  containing  a  quasi-rigid  rotor  and 
arigidplatfbnn.  Itisobserveddiatdieadectioooftfaebody  tobetherotoranddiebody  to 
be  die  platform  is  completely  srintraiy ,  and  no  physical  distinction  exists  boween  d^  two 
bodies.  Therefore,  by  analogy,  die  equation  describing  the  torque  on  die  quasi-rigid  rotor 
as  a  function  of  rotor  vaiiaUes  inust  be 

Tmt  -  <«) 

Now  let  Tpiy  and  Ti/^  represent  the  motor  torques  on  the  platform  and  on  the  rouv 
respectively  for  the  system  containing  bodi  a  quasi-rigid  rotor  and  a  quasi-rigid  platfonn. 
Then  the  separation  axiom  would  require  die  fcdlowing  two  conditions 


TplM  =  ^ 

X 

(67) 

Tr/m  =  Trim  «  ^+/,ov 

A 

(68) 

Once  again,  siiKC  die  system  is  an  actkm  -  reaction  pair 

TpiM"^TpiM  *  (^ 

(69) 

andthra 

-  +  —*-(/, 6tt’¥la  »/)  B  |li|  Osin  9 

A  X 

(70) 
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TteiMbiHQroilBfiondkaHestfiattf  itaMinxaDorbeii^ttive,tluis 


^£+iiso  .) 

A  A' 


This  is  seen  to  be  die  existing  eneigy  sink  criterion.  Therefore,  despite  the  presence  of 
motor  torque,  the  existing  stst^ty  criterion  sdU  lilies.  Referring  to  Equntions  (67)  and 
(68),  die  second  term  in  die  equations  would  rqxesent  the  torque  ^^lied  by  die  damping 
tnafthMii«in  to  die  platform  and  the  rotor  reflectively.  Equations  (67)  and  (68)  state  that  the 
motor  torque  minus  the  damper  torque  will  equal  the  net  torque  of  the  platform  and  the 
rotor  respectively. 

In  determining  the  revised  energy-sink  stability  theoiy,  the  energy  dissipation 
equation  for  die  system  with  no  energy  contribution  from  the  motor.  Equatim  (37),  must 
be  rewritten  to  account  for  the  motcv  torque.  Thus 


^  total  +  1^  =  —lfXd>s-‘l$*X'  6>s' 

Substituting  in  Equations  (67)  and  (68),  one  arrives  at 


^  *  ^Dtotai*^  “ 

=  6d  +  to  -U  ^p/m) 

By  assigning  the  motor  torque  values  as 


*  Ta/iir  =  -’^p/M 

dien  die  above  equation  can  be  rewritten  as 

£  '  iDtotal  +  ^  «  ^d*^d*Tm(X'-X) 


01) 

(72) 

(73) 

C74) 


Because  die  total  time  rate  of  change  of  energy  dissipation  equals  the  rate  of  change  of 
energy  dissipation  in  the  rotor  plus  diat  in  the  ftiatfmm,  die  rate  of  wmk  due  to  the  motor 
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tocquemustbe 


W  m  •TmO  OS) 

Referring  to  the  example  from  Chiq)ier  1,  the  prolate  dual-spin  communications 
satellite,  typically  the  rotor  nutation  fiequency  would  be  n^ative  and  the  plalfonn  nutation 
fiequency  would  be  positive.  For  rotor  spin-up  Ty>Q  and  the  motor  torque  would  add 
energy  to  the  system,  and  for  rotor  spin-down  Tm<0  and  the  motor  torque  would  remove 
energy  from  die  system.  For  the  arbitnuy  system,  die  motor  torque,  Ty,  and  the  sign  of 
die  term  A'  -  A,  would  determine  whedier  the  motor  adds  or  removes  energy  from 
the  system. 

The  rates  of  change  of  the  energies  of  die  system  may  now  be  represented.  Rewriting 
Equadon  (46)  to  determine  die  total  energy  dissipation  of  die  system 


taul  —  i 


06) 


The  rate  of  work  due  to  the  motor  uxque  can  be  expressed  by  combining  Equations  (75) 
and  (68)  (or  Equation  (67) )  to  arrive  at 


X 


(77) 


The  rate  (rf  change  of  the  kinetic  energy  of  the  system  as  a  function  of  platform  core  energy 
and  system  parameters  was  determined  previously  as 

£  =  £c'  +  h  +  o)  +  /j  0)5'  a  (52) 

Substituting  the  above  two  equations  into  Equation  (76),  the  expressim  for  total  energy 
dissipation  of  the  system  may  now  be  written  as 


total  ~  £c'  +  h  <f{^s  +  ^  0>s  O’-—  O-/,  <05  O 

X 


OS) 
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BacttMe  0^' 4- a«  dds  cqnidQn  cn  be  Mier  sfaBiriified 

totcud  -  (79) 

GoaqMiing  this  eqiiatioa,  rqnesenting  the  system  widi  a  quasi-rigid  loior  and  a  quasi-rigid 

platfonn,  with  Equation  (S4),  representing  die  system  widi  a  rigid  rotor  and  a  quasi-rigid 

* 

platform,  reveals  an  additional  tetm,-^  g  The  first  term  of  Equation  (79)  represents  the 

X 

rate  of  change  of  core  energy,  with  die  platform  as  the  core  body.  The  second  term  will 
account  for  die  change  in  energy  associated  with  a  change  in  the  relative  rotation  rate  of  the 
rotor  with  respect  to  the  platform.  The  final  term  accounts  for  die  energy  loss  due  to  the 
quari-rigidity  of  die  rotor  that  is  not  represented  in  die  platform  core  body  expression. 

The  case  that  will  be  analyzed  is  that  which  occurs  when  diere  is  a  constant  reladve 
rotation  rate  between  the  platform  and  the  rotor.  For  die  remainder  of  the  analysis,  let 

^  =  0  (80) 

and  the  total  rate  of  eneigy  dissipadon  becomes 

=  6c'-—  (81) 

X 

Further  simplification  can  be  achieved  by  noting  that  as  X'-X  and  fiom  Equation  (37) 

Therefore 


£i)  +  4i)'  *  — “j 


which  reduces  to 


+  (83) 

X  X'  X' 

A  similar  development  can  be  performed  using  rotor  core  energy  vice  platfmm  core  energy. 
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Theiesukis 


m  ^ 

X  A'  *  A 


(84) 


Tliere  is  an  expected  symmetry  between  Equations  (83)  and  (84).  due  to  the  arbitrary 
assgnmentofonebody  of  die  system  as  die  rotor,  and  die  odier  body  as  die  platfonn.  To 
confirm  the  results  of  Equations  (83)  and  (84).  one  must  prove  that  udieno  isconstant 


A  A' 

Ihldng  the  derivative  of  die  rotor  core  energy.  Equation  (48) 

&C  *  f|  tcuU  rOi  +  ta$at 

*  f|  Mo/  4*1  h  tauli^s  +  +  ci) 

and  sitnilatly,  for  the  platfonn  core  energy 

local  4*1  +  f#  locol  4*/  4*/ 

Substituting  Equation  (87)  into  Equation  (8^  and  noting  that  0’s  0. 

+  Is  tetai<f^s 

Multiplying  dnou^  1^  the  platfonn  nutation  fiequency 

&C  A'  «  ic  A'  +  /,  total  A'  0®/ 

and  recalling  that  A' s  A  0, 

£cA'  *  A  +£c'O’+/iiocalA'0®/ 
s^ndi  is  die  same  expression  as  Equation  (8S)  if  it  can  be  proven  that 

ic  +Is  total  ^s'  *  0 

EHmiMting  die  trsnsveise  angular  velocity  of  Equation  (47)  by  substituting  in  Equation 


(85) 


(86) 


(87) 


(88) 


(89) 


(90) 


(91) 
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(3Q.<jaeaBivefat 


,2 


Taidiig  the  derivative 

tc’  -  .fe <«‘.t '»>»'), <»•  * >•' o; imm 

h  talat  h  toui 

and  rimplifying 


(92) 


(93) 


m  Ois  dtf* Igig  Wg  6i<  fg  Oig  dg  )  ^  /j  total h  total  fOs 

A  tMol 

Noting  that  Igiattd  •  Is  +  h'  and  recalling  diat  dv  «  ov'  because  «  0,  the  above 
equation  can  be  reduced  to 

Ec  *  total  .  -1 1  ii"..  (V5) 

•t  tetai 

finaUy,  invoking  die  definition  of  the  idatfonn  nutation  fiequency.  (Hie  airives  at 


U  " -UmX  i>,  (96) 

Hieiefore,  Equatioa  (91)  bolds  and  die  revised  staUlity  oiierion  can  be  expitssed  as 

^  ~  so  (97) 

A  A  y 

A  few  remaiks  can  be  made  ctmceming  this  stalnli^  criteri(Hi.  The  third  expressicm  is  the 
existing  energy  sink  stability  criterkm.  This  critericm  must  equal  die  stability  criteritxi  as  a 
function  of  the  rotor  core  energy  whi<di  must  equal  die  s&dxliQr  criterion  as  a  fiinctxHi  of  the 
platfonn  one  energy.  It  can  be  seen  that  one  no  Icmger  needs  to  know  the  energy 
disripatiiHi  rates  in  die  platfonn  and  in  the  rotor  to  determine  stability.  By  knowing  or 
postulating  the  cent  energy  over  time,  the  stability  of  the  system  can  be  deteimined. 
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CoBtinuing  with  die  piolMte  dual-spinner  example  of  Chapter  1,  any  one  of  the  above 
eiqiressioiis  apply.  Fbr  a  staUe  system,  die  rotor  core  eneigy  wiU  be  positive,  and  increase 
over  time.  Additionally,  the  rotor  natation  frequency  will  be  negative,  resulting  in  a 
negative  eiqiressioa  for  the  rotor  core  energy  stability  criterion.  Hie  platform  core  energy 
will  be  positive,  but  will  decrease  over  time.  The  platform  nutation  frequency  will  be 
positive,  resulting  in  a  negative  ex|»ession  for  die  platform  stabili^  criterion.  Acocrdiiig  to 
Equation  (97),  die  platform  and  rotor  staUlity  criteria  must  equal  one  anodier.  From  the 
numerical  simulation  of  a  dual-^nn  quasi-rigid  axisymmetric  system,  die  rotor  and  platform 
core  energies  as  a  function  of  time  will  be  determined  and  graphed. 

D.  NUTATIONAL  MOTION 

The  developroent  of  a  modified  expression  for  the  nutation  angle  as  a  function  of  time 
may  now  be  presented.  The  actual  nutation  angle  of  the  system  is  defined  as  0.  The 
nutation  angle  as  a  function  of  platform  core  energy  will  be  represented  by  rji',  and  the 
nutatkxi  angle  as  a  funetkm  of  rotor  core  energy  will  be  represented  l^r).  The  derivation  is 
similar  to  the  one  prev^sly  given  in  this  chiqpter,  except  that  the  total  energy  of  the  system 
has  been  replaced  by  the  core  energy  of  the  system  to  eliminate  the  transverse  angular 
velocity,0)^.  The  derivation  for  the  platform  core  energy  will  be  shown.  From  before,  the 
angular  momentum  of  the  dual-spin,  quasi-rigid,  axisymmetric  system  is 

*  A  w*a/ +  (a  0>s  +  //  a/)*  »  0>i^ +  {ls  <r+ h total  0>/Y  (26) 

Cjombining  this  with  the  platform  core  energy,  and  solving  for  the  platform  angular  velocity 
about  the  spin  axis 

ft  teul{ls  U)tot'~  h  lotad^*  —  2 io*at  A  +  A*  +  2  £c^ /» ima/ —  h*  ■  0  (98) 
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O'  -  (loo) 

\  hiotatl  •tmai'  Mttquii 


Tile  imtialooaditioiis  at  r»OwiUdeiBniiiiie  die  proper  sign  of  Equatkn  (99).  As  before, 
ooMiiiiiity  win  ensure  dussigD  for  all  t>0.  Substituting  Equation  (99)  into  Equation  (11) 


rt  .  .  CO.-.  ao.) 


The  timexaie  of  change  of  die  nutation  angle  wiU  detennine  the  sttbility.  Differentiating  the 
above  equation  results  in 


if  >  1 

si«n  *  (±v^) 

A  stable  system  would  require  that  die  lower  sign  be  chosen  for  the  radical,  thus 

When  this  sign  selection  is  ^iplied  to  Equatkm  (S2),  one  gets 

v>  toco/' 


(102) 


(103) 


(104) 


This  can  be  rewritten  as  die  well  established  dual-spin  stability  condition,  as  written  by  P. 
C  Hu^es  [Ref  6], 


{ft  total  ■"  1$  local)  '^igO  ^0  (lOS) 

It  is  impoftant  to  note  duu  this  aiudysis  does  not  produce  a  contradiction  to  Equation  (41). 
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r.  die  nxMified  ouottkxi  aa^  OB  be  derived  with  leqiect  10  the  rotor  at 


Theeefoie.  if  ooe  is  givcfi  die  core  eaergy  or  die  potftilatnd  one  energy  as  a  fiinctioa  ci 
time,  die  notadoiial  motiao  can  dien  be  deiennined.  Tins  leads  to  an  extremely  iinpoftant 
conclusion.  By  deiennining  a  suffideady  accurate  postulation  of  die  core  energy  of  a  dual- 
qiin,  quasi-rigid.  axisyimnetric  spacecraft  over  time,  ooe  can  predict  die  nutational 
behavior  and  die  stability  of  dutt  qMcectaft.  Numerical  simulation  of  such  a  system  will 
confirm  dus  conclusion. 


m.  DYNAMICAL  EQUATIONS 


A.  MINGORl  DUAL-SPIN  SYSTEM 

AqwdficoioddisieqiiivBdtovalidMtheiiiodifialeDafSMniAMbilttythe^  Hk 
devdopoett  of  the  pitvkws  chiiicen  was  com|rtetely  fenenl  in  nature  and  q)plies  to  any 
dnal-qnn  q»oecnft  widi  an  arbitraiy  ffanyinf  mechaniam.  D.  L.  Mingori's  dual-qra 
syston  [Ref.  3]  profvided  the  needed  model  required  to  validate  die  proposed  stability 
theory.  Figure  (3).  Additionally,  the  conoplete  non-linear  equations  of  modon  were 
presented  by  Mingoii;  however,  the  expressions  for  the  important  parameters  in  attitude 
dynamics,  angular  momentum  and  kinetic  and  total  energy,  were  not  presented  in  his  pqier 
and  had  to  be  derived  before  any  analysis  could  be  performed. 

The  Mingori  S3rstem  is  comprised  <rf  two  symmetric  rigid  bodies,  the  lower  which 
shall  be  defined  as  die  rotor,  and  the  upper  botty  shall  be  defined  as  the  platfonn.  By 
convention,  all  terms  refoning  to  the  platfoim  will  be  die  same  notation  as  that  of  the  rotor, 
excqn  that  dieywiU  cany  die  prime  maik.  Both  the  rotor  and  die  platform  have  coordinate 
axes  fixed  to  die  body  and  located  at  die  leqiiective  centers  of  mass.  The  distance  between 
die  centers  of  mass  is  qiecified  by  L.  The  entire  q^aceciaft  has  coordinate  axes  fixed  to 
the  qiacecraft  center  of  mass,  denoted  by  die  douUe  prime  coordinate  axes,  and  rotating  in 
die  same  manner  as  die  rotor  coordinate  axes.  The  coordinate  axes  bj.  and  b3''areall 

Gcdlinear.  The  qmcecraft  center  of  mass  will  vaiy  along  the  bs  axis  as  the  point  masses  in 
the  rotor  and  the  platform  oscillate.  The  only  relative  modon  die  platform  widi  reflect  to 
die  rotor  is  angular  rotation  about  die  ba  axis.  A  motm’  driven  by  a  control  system 
maintains  a  constant  relative  rotation  rate  o.  The  angle  between  a  line  parallel  to  bi  and 


passing  through  the  platform  center  (rf  mass,  and  bi' is  rqnesenied  by  yr »  at. 
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Ikfii^cn  Daal-Spin  Qiiasi>Rigki  Axi^naneok  Sysi^ 
Figure  (3) 
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PlliJUm  .  :,JI, 


Banek  body  containt  a  mass-qning-daahpot  mechanism.  An  actual  spacecraft 
onderfoes  damping  from  various  mechanisms.  Undesiied  damping  can  occur  due  to 
structural  elastid^  and  liquid  propellant  slo^.  Nutation  dampers  are  incorporated  into 
spacecraft  to  inqnove  the  nutational  motkm.  The  dual-spin  axisymmetric  system  with 
mass-spfing-dasIqxM  dampers  cannot  accurately  model  an  actual  spacecraft,  but  is  used  to 
illustrate  and  validate  die  dieoiy.  A  description  of  the  rotor  is  as  follows.  The  mass- 
qpring-dashpot  mechanism  can  be  modelled  by  a  particle  mass  mi  attached  to  a  firing  with 
constant  k  inside  a  tube  filled  with  viscous  fluid  widi  damping  coefficient  c.  The  moticMi 
oi  the  particle  is  constrained  parallel  to  the  ba  axis  only.  At  rest,  the  particle  mass  lies 
along  the  rotor’s  bi  coordinate  axis,  at  a  distance  a  from  die  rotor’s  center  of  mass.  Three 
balancing  masses,  m2,  m3,  and  1/14,  each  equal  to  the  mass  of  the  mass-qning-dashpot 
mechanism,  are  ripdly  fixed  a  distance  a  on  tl«  b2,  -bi,  and  -b2  axes.  These  masses 
maintain  the  axisymmetry  of  the  system  about  die  bs  axis.  Displacement  of  the  particle 
mass  mi  is  denoted  by  the  variable  2.  A  siaqilification  in  this  paper  of  the  Mingoii  dual¬ 
spinner  system  is  the  assunqidon  of  a  massless  spring-dashpot  system.  Thus  the  particle 
of  the  mass-spring-dashpot  system,  mi,  is  the  same  mass  as  the  correqxniding  three  other 
balancing  masses,  m2,  m3,  and  1x4.  The  platform  can  be  described  in  a  similar  manner, 
widi  all  notation  modified  with  the  prime  mark. 

The  dual-qiin  system  center  of  mass  ccxxxlinate  axes  rotates  at  the  sanoe  rate  as  the 
rotor  coordinate  axes,  but  is  located  akmg  die  b3  axis  at  a  distance  Zcm  from  the  rotor  center 
(tfmass.  This  distance  will  vary  as  the  particle  masses  are  di^laced.  Relating  the  dual- 
^lin  system’s  coordinate  axes  to  the  rotor  coordinate  axes  was  arbitrary.  Equivalent  results 
would  be  adiieved  by  selecting  die  platform  coordinate  axes  instead. 
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B.  ANGULAR  MOMENTUM 


The  anpdar  mcmieittiim  of  the  Ikffinfori  dnal-^piiu  qoiSMigid,  axisyina]^^ 
nofw  deiived  fian  film  priflcqiles. 

The  angnlar  momentum  (moment  of  momentum)  of  a  panicle  of  mass  m  j  located  in 
botfyiB  is  defined  as 

hi-riXmiV  (1(») 


where  hi  is  die  angular  momentum  of  the  i  tfa  pardde  widi  respect  to  die  system  center  of 
mass,  riqiecifies  the  location  of  the  ith  particle  with  reflect  to  the  system  center  of  mass, 
and  V  »  is  die  abscdute  velocity  of  die  i  di  panicle  widi  reflect  to  the  inenial 
reference  fiameM  Expressing  the  absohite  vdoci^  in  the  system  reference  firame 


N;*eM 

V  +ri  + 


N-^B 

at  X  ri 


(109) 


adiere  B  represents  the  reference  coordinate  axes  of  the  system.  For  die  following 
derivation,  all  displacements  and  velocities  are  referenced  widi  reflect  to  the  rotor  (hi) 
coordinate  axes.  The  angular  momentum  vecmr  is  rewritten  as 


iN^Sem 

hisrjxmi\  V 


N^B 
+  fi  +  ®  X  r 


.1 


(110) 


^iplying  dus  equation  to  particle  1  widi  mass  mi  andposidon 

n  bi +  0  b2  +  (z-:^)b3  (111) 
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0  -a{2-Xa^  »1 

0  +  0  ^  * 

-fl(*-*aii)  0  tf*  \\^ 

0  -{z-xgm)  0  1/  v«,,  \  r  0 

(*-*«)  0  -a  Vo,,  }+  -a(i-ieJ 

0  a  0  J\  Vamt  I  L  0 

Rr  the  rigid  body.  Equation  (1 10)  is  ^ipfied  to  1  differential  volume  at  location 

r^*u  bi +  v  b2  +  (w-2c/ii)b3  (113) 


Integrating,  one  anives  at 


/i  +M  zii  0 

0  h 

0  0 

0  MZcm  0 

-MZan  0  0 

0  0  0 


(114) 
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SigailarftalftwitWw  m  perCnmied  ftir  the  tematiiing  tevCO  point  mailBS  ind  ifae  plltfonn. 
For  the  RNor.  one  Brives  at 
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(116) 


hie 

/ 


(/-!«)*♦ 


(2tf^  +  4(/-*«)f‘^ 


-ii^a'i*Go$(<rt)  a"  t'm{at) 


-  «i'a'x'oof(9r) 


/s' ♦4  m' o'* 


«! 

«2 

103 


\ 


♦ 


+ 


-  m' o' <r  cos((r  f)  s' ♦  in' o' sin(9 1)  s' 

-  in' o' <r  sin(a  f)  s' -  m' o' cos(9 1)  s' 

(/3'  +  4m'fl^)ff 


0 


-(4/'  +  4m')(/-s„) 
-in's' 


(M'  +  4m')(/-Sc.) 

♦  m's' 


0 

0 


0 


0 


0 


y 


/ 


The  angular  momentum  equations  have  terms  conesponding  to  the  velocity  of  the 
center  of  mass.  Angular  momentum,  when  taken  about  the  center  of  mass  of  a  system,  by 
definiticm  must  be  indqiendent  of  the  translation  of  the  system's  center  of  mass  (with 
leqaect  man  inertial  reference  fiame).  To  verify  that  tins  occurs  in  die  above  equations,  the 
equation  for  the  position  of  die  center  of  mass  is  substituted  inm  the  center  of  mass  velocity 
tenns,  and  then  these  equations  are  equated.  If  they  have  die  same  magnitude  but  opposite 
sign,  diey  will  cancel  each  other  out  and  it  will  be  proven  that  the  system  angular 
momentum  is  indqiendentof  die  trandation  of  die  system  center  of  ma^ 
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One  mast  prove 


Looldiig  at  die  bi  oooqmnents,  one  finds 

-[(iW'  +  4  m')(-  /  +  Son)-  m'  /]  *  [(Jl#  *4m)zam-fnz] 
Rewriting  dus,  one  obtains 

{il#'  +  4m')/  +m2  +  m'i'*(Af +  4m  +  Ar  +  4m')2i», 
ThecenierofmasswidiieferenoetDdiebi.  b2.  bs  ooocdinate  axes  is 

2  Mi  Zi  +  mi  Zi 


(117) 


(118) 


(119) 


(120) 


mz  +  (M'  + 4  m')/  +m'z* 

Ml  4*  iR|  +  4  HI  ^M^  +  4  HI 

i 

Substituting  Equation  (120)  into  Equation  (119)  results  in 

{M'  +  4m')/  +«z  +  «'z'  I  (121) 

\  Jlf  +  4m+M'  +  4m'  / 


This  reduces  to 


(at  +  4 m') /  +  m z  +  m' z'  «  (aT  +  4 m') /  +mz  +  m' /  (122) 

Therefore,  die  angular  momentum  along  the  bi  axis  is  independent  of  the  velocity  of  the 
center  (tf  mass.  Verifying  this  along  the  other  axes  can  be  done  in  a  similar  manner.  This 
proves  that  the  angular  momentum  of  the  entire  system  about  the  sysmm  center  of  mass  is 
indqiendait  of  the  oanslaticm  of  the  system  center  of  mass.  The  angular  mmnentum  of  the 
entire  system  can  be  written  as 

*  l»j#+4ii»+ iijr+4ji<'  (123) 
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C.  ENERGY 

The  kinetic  energy,  £,  and  the  total  energy,  Eunai*  are  now  derived.  From  first 
piincq)les  die  Idhedc  energy  of  a  difierendal  particle  is  written  as 

(126) 

where  the  absolute  velod^  is  defined  fiom  before 

Nriim  N^cm  N-^M 

V  »  V  0)  Xfrf*  (109) 

The  Idnedc  energy  of  die  Mingoii  system  is  determined  by  integrating  die  differential 
kinetic  energy  over  the  rotor  and  over  the  platform,  and  summing  the  differential  kinetic 
energy  equaticm  over  the  eight  pmnt  masses  to  arrive  at 

r 

i  L  ( ( ''^ )  •  ( ''v*' )  1  <*"'  + ^  1 1 V )  •  ( V ) ) 

Pdfocming  the  steps  on  particle  mi,  the  following  is  obtained 

(128) 

A  V’)+{ri  ^i)+({''5"’x'’>)- (''S“x'’i))+  1 

^  i- 1  1 2  •  ri)+  2  ("^  •  ("S'-x  r,) )+  2  (r,  .  ("^'-x  r,) )  j 

Substituting  in  the  values,  tme  arrives  at 

z  +  Vjn  y  +  *  +  (i  “  2cmF  +  ~ 

a^a}}  +  (z-ZamfiOi-2a  (r - Ze*) Oh  013  + 

a^O>i  +  2(zVemz-iemVcmt)+  (129) 

2lz(02Vemx-Zcm0>2yemx+aa>3Va„y-\ 

1  *<»1  Vcmy  +  Zema>lVemy-aO>2ycmi  I 

2(-aza)2  +  azcmCt>2) 
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loll  oponnoB  n  nn  pmniiisa  ow  aw  oowr  noac  puigci  in  aw  nnor.  mranaaig 
tlww  now  siqis  oa  tfw  boify  of  the  nNor. 

(130) 


Eli' 


i|_(2r 


V***  •  ^^)+(ri  •  ri)+  { (^3* X  n ) .  (^2«» X  n) ) + 

•  I‘i)+ 2  (^V^  .  r,) )+ 2  (h  •  (^3‘*X  n) ) 


dM 


Substituting  in  values,  and  qiediying  tfw  position  vector  of  the  differential  particle  as 

r^«iibi  +  vb2  +  (»v-*oii)b3  (131) 


then 


(132) 


J  Mmr 


yimx-^yLy+^x  +  iim+**'^0>2-2wZcm0)i  + 
ximO)2+V^ioi-2wV  i»2^  +  2vZcmO>2a>3+W^a)f  + 
Z^i0i+U^i03’-2w  2cmi»i-2uw0i  (»3  +  2U2cm('^l 
V^tuf  >2  Vtl  0)1  0)2 -i'U^  0)2**’ 2 (-Zcm  Vcmt)  + 
2(*^<»2Vo«,-ZciiiO)2Vo,,-VO)3V«.,-WO)iVa„,+\ 

\  ZcmiOlVcm,  +  UO)3Vc„y  +  V(»lVamt-ua>2Vemt  I 

2(-V2em0fi  +  U2em  tOz) 


Integrating  over  die  limits  of  the  symmetrical  body,  the  kinetic  energy  of  the  rotor  becomes 


Em-^U 


Vi,  »  +  Vi, ,  +  vj, , + i|,  +  li,  (<»f + a>|) +' 

/l  0)f +  72  0)2  ■*’/3  0)3 - 
L  2  2cm  0)2  Vcmx  +  2  Zcm  0)i  Vcm  j~2  Zem  l^cm  x 


(133) 


These  same  steps  are  performed  on  the  platform  and  rm  the  remaining  seven  particles. 
These  equations  are  dien  ctxnbined  and  simplified. 
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Tile  my  Idoedc  eacfgy  to  die  Miofari  qracm  is 


^M««,(v2,,+  vL,+v2,,)+ 
l(/i  mi+l2€»i+h  »f  *l2ai*  h'ian  +  o)?) + 

|■(z^ { Aif  oij)  -  2  a  z  Ah  013) -f 


(134) 

\ 


^(z'^ -fz^ (a>i  oil) -let ^ oos(<Tr) a>i  a>3 *>' 2 o' z' siii(ar) 012  ob) '<* 

m{-az€^* 

nt  (-et  coi{at)  loi'^eti*  assert)  t»i-et/  (Tcot(at)  Oh-iti'  asu^fft)  02} 

1  («>f  +  <02)(^  z  +  m'  z'f  -{m  z  +  m'  z'f 
*2  Mt 

^(<»f  +  AlfK^ ^ ^ (^Mj  ”F ” 

U2  2\iin2*iit  /  *{hf  •¥4m')l\ 

- Wr - -}- 


\ 


nt  2* 


The  total  energy  of  the  Mingori  system  can  be  determined  by  adding  the  kinetic 
eneigy  of  the  syston  to  die  potential  energy  of  the  system.  The  only  potential  energy  of  the 
system  is  the  energy  stored  in  the  springs  of  the  mass-spring-dashpot  damper.  From 
Hooke's  law,  die  system  potential  energy  is  written  as 

(135) 


The  total  energy  of  die  h^ngori  system  is  then 


D.  M1I^KKHU*S  EQUATIONS  OF  MOTION 

li(6il|ori's  eqoatioiis  of  moikn  [Ref  3]  are  wrinen  as 


-2  (C+ -  ®2  «3) + i  Wi]  + 
z  [-2  C  (Oi-f  2  i  o»i  +  z  (<bi  -  0)2  o^)-o(<^  *  <>>10)2)] 

-  2  (C- /l)k -  <»2«3)  +  »ll  + 

z'[-  2  +  2  0)i  +  z'Clbi  -  0)20)3)- o' ootvr(<^  +  0)10)2)]  + 

o'sinvr  {z'  +  z'[{ob  +  of-o^]| 

A  0^— (C— A) 0)i  0)3  —Ji  ooni  +  2MT(f^oi^-¥Mx C^{ci>2  +  0)i  0)3)  + 


m 


-2(C+l^[z{d>i  +  o>i  0)3 
Z  [-  2  C  0)2+  2  z  0)2  +  z  (<»2  +  0)i  0)3) 


+  Z0)2]  + 

-a[z  +  z(o)J-0)?)]/^ 


(137) 
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-  2(C-  /i)[z'(«2  +  0)10)3)+ z'  0)2]  + 
z' [- 2  C  <»2  +  2  z' 0)2  +  z' (o)2  +  <»i<»3)  ] - 
\  o' oosy^fz'  +  z'[  (o)3  +  -  o)f] } -o'  smyrz'(a)3  -  o)i  0)2) 

C  o^  -  m  o  [2  z  o)i  +  z  (a)i  -  0)2  0)3)] - 
m'o'  sinyr[2z'o)2  +  z'(o)2'*'<Vi  0)3)]- 
m'o'  cosyf[2z'o)i  +  z' (01-0)20)3)]  »0 

m(l  -p)z-m'pz'-mo((^-0)i  0)3)- 
m  (o)i  +  0)1 )  [z  ( 1  -  p)  -  ^  -  p' z']  +  c  z  +  *  z  »  0 


-m  p'z  +  m'(l  -p')z'  + 

m'o'(sinyr[o)i  +  0)2  (0)3  +  2  o)]- cos  yr[c^- 0)1  (0)3  +  20)]  )- 
m'(o^  +  o)|)[z'(l  -p')  +  /i  -pz]  +  c'z'  +  it'  s' «  0 
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CSiapier  IV  disoisses  how  tfiese  equatioiis  are  adapted  for  nse  in  the  numerical  integration 
routine. 


E.  NUTATION  ANGLE 

The  nutation  angle  of  die  dual-qiin  system  can  be  determined  by  substituting 
Equations  (124)  and  (12S)  into  Equation  (34)  to  anive  at 


d*cos“* 


COS“* 


(139) 

(-m  a  2  -  m'  a'  s'  co5(<rt))  a>i  +' 
(-m'  o'  z'  sin(<Tr) )  <»2  + 

{13  +  4  m  +  73'  +  4  m'  o'*)  013 
+  (73'  +  4m'o'*)<T 
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F.  CORE  ENERGY 

The  {riatfonn  core  energy  and  the  rocor  core  energy  were  defined  in  Equatioos  (47) 
and(48).  The  nutation  angle  as  a  fiinctioa  of  core  energy  was  then  developed  and  is  stated 
in  Equations  (101)  and  (106).  These  equadon  can  be  used  to  predict  die  nutation  angle  of 
die  system  as  a  function  of  time.  Because  the  core  energy  as  a  function  of  time  is  not 
avaUaUe  for  the  predictioii,  a  postulated  core  energy  must  be  developed.  The  initial  value 
the  actual  core  energy  and  the  postulated  core  energy  must  match  and  is  dictated  by  the 
inidal  conditions.  Fsnmeters  (tf  die  postulated  core  energy  must  be  selected  to  accurately 
model  the  actual  core  energy,  to  include  the  final  energy  state  and  the  rate  at  which  the  it 
approaches  the  final  energy  state.  Two  models  were  considered. 

1.  Exponential  Core  Energy  Model 

An  expmiential  representadon  of  rotor  and  platform  core  energy  as  a  function  of 
time  are  expressed  as  fdlows 

£c  yoiriifafrf  W  *  (Eco  -  Ec final)  eh'^)  +  Ec final  (140) 

Ec poanlated  W  *■  (Eco' - Ec final) +  Eq final  (141) 

The  initial  core  energies,  E(^  Eco't  are  determined  by  the  initial  conditims  of  the  system. 
The  final  core  energies,  Ec  Au/,  Ec  finals  as  weU  as  the  exponential  factors  r,  /,  must  be 
selected.  The  medioddogy  for  selecting  these  values  is  explained  in  Chtqiter  V,  Con^uter 
Analysis. 

2.  Verhulst  Logistic  Core  Energy  Model 

The  exponential  model,  as  will  be  explained  in  Chapter  V,  has  excellent 
agreemoit  for  stable  cmditimis,  but  performs  poorly  for  the  unstable  conditions.  As  an 
altenutive  model,  die  following  first  order  dififerential  model  is  selected  for  the  rotex’  and 
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^  pcmaltiadi^ 


rw 


. _ _ 

Ect  +  \Ec/kml-Ec^^"'*^ 

Ec:AEcfi^'-BcM-^*^ 


(142) 


(143) 


This  type  of  equatioa  was  fim  introduced  by  P.  F.  Veihulst  lo  model  human  and 
odwrpopulatioas  [Ref  7].  It  is  often  refentd  to  as  the  Vesfaulst  equation  or  the  logistic 
equation.  AMxxii^  population  dynamics  tqtpears  to  be  unrelated  to  the  stability  ttf  a  dual- 
qpin  system,  the  behavior  over  time  ci  this  equation  conqiared  to  the  stable  and  unstid>le 
systems  provides  some  insighL  Hgure  (4)  shows  die  Veriiulst  equation  versus  time  for 
varying  initial  conditions. 


Rgure(4)  Veihiilst  Logistic  Equation  Versus  Hme 
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It  can  be  noted  diat  if  die  initial  value  of  the  variable  0n  our  case  core  energy)  is 
jreaterdian  the  equiKhrium  value,  it  will  approach  die  equilibrium  value  in  a  manner  sinular 
to  diat  of  an  exponential  decay.  Fbr  initial  values  that  are  less  than  but  within  one  half  of 
the  equilibrium  value,  die  variable  will  asynqittitically  ^iproach  die  equilibrium  value.  If 
the  initial  value  is  less  than  one  half  of  the  initial  value,  die  variable  over  time  has  a 
somewhat  difieient  shape  as  it  appioadies  equilibrium.  Initially  the  skqie  is  very  small,  but 
increases  to  a  maximum  at  about  the  one-half  of  die  equilibrium  position.  After  this 
inflection  point  die  vaiiaUe  qiproaches  equilibrium  asymptotically.  The  value  of  Nq  equals 
if  is  an  asymptotically  stable  equilibrium.  A  value  of  No  equals  zero  is  an  unstable 
equilibrium.  If  the  value  of  No  is  slightly  greater  than  zero,  then  Nit)  will,  as  t  -><>•. 
achieve  the  staUe  equilibrium  of  K.  These  types  trf  curves  will  be  utilized  to  describe  both 
the  stable  and  unstable  dual-spin  system.  The  initial  core  energies,  Ecg.  would  be 
determined  by  the  initial  conditions  of  the  system.  The  final  core  energies,  Ec  Aia/> 
£c  final*  must  be  known  or  a  best  estimate  used.  The  exponential  factors,  r,  /,  must  be 
determined  experimentally,  and  are  a  function  of  the  system  parameters  and  initial 
conditions.  Chiq)ter  V  provides  additional  details  on  die  selection  process. 


IV.  NUMERICAL  SIMULATION 


A.  NUMERICAL  SIMULATION  EQUATIONS  OF  MOTION 

Tlie  equadons  for  the  dynamical  quantities  needed  for  analysis  of  the  dual-spin 
system  were  previoosly  developed.  Man^polatioo  was  leqniied  to  make  these  equatkms 
suitable  for  numerical  analysis,  and  is  described  bdow.  Tlieoompuier  code  is  included  as 
^ipemfixB. 

1.  Minfori's  Equations  of  Motion 

The  equations  of  motion  for  die  dual-q»n  system  are  listed  as  Equation  set  (137) 
and  (138).  The  Runge-Kutta  numerical  integration  routiiie  necessitated  that  the  equations 
of  motion  be  a  set  of  first  Older  differential  equations.  Madiematical  manipulation  was 
lequiied  to  imt  diem  in  diis  foim.  VaiiaUes  iqnesentiqg  groups  ai  ttxtas  (A^  Bh  Ch  etc) 
were  introduced  to  sinqilify  the  manipulations  of  die  equadons  and  provide  a  suitable 
format  for  iiKxxpoiating  the  equadtms  into  die  conqioierprogriun.  Mingori's  equadons  can 
be  eiqiiessed  in  terms  of  these  variables  and  the  five  dme-dqiendent  variables  of  motion 

^26  ’►  ^27  +  ^24  +  ^28  *  0 

B25  ®2  +  F23  +  Ri3  Z  +  B2\  +  B26  *  0 

Ciomi'*'Cs(^'»-Cc^4-Cn«0  (144) 

Di  z + Z>2  +  f>8  =  0 

£i  z  +  £2  J +  £5  ffli  +  £io  «  0 

where  the  variables  Ai,  Bi,  Q,  Di,  Ei,  Fj,  and  Z,-  are  defined  in  the  Notation  section. 
Qeariy,  these  equations  are  highly  coupled.  Throu^  a  series  of  manipulations,  including 
substitution  and  combining  sets  of  equations  to  ehminate  common  variables,  one  can  arrive 
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l-^10^2I,gg.ga 

C  A26  C  B25 


A26  ^  A26  A26 

^25  ^  ^25  ^25  ^25 


The  equations  remain  coupled,  so  the  sequence  in  which  the  equations  are  numerically 
int^rated  is  important  Because  die  equations  for  z  and  the  equation  for  z' are  functions  of 
Zands',  z  and  z' must  be  integrated  fost  Similarly,  z  and  z' must  be  integrated  before  0)3, 

013  before  tUi,  and  tui  before  013.  In  the  conqiuter  program,  the  seven  variables  that 
describe  die  motion  are  stored  in  a  seven  column  matrix  wliere  a>i »  y[i][j],  02  «  y[2]  [ il, 

03  -  yI3]Q],  z  «  yI4]Ql,  z  »  ylSlDl.  z' «  yieiQ].  and  s'  -  y[7]Dl.  The  index  i  identifies  the 
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matrix  location  for  each  set  saved  vaxiabtes  over  time.  All  other  needed  quantities  can  be 
calculated  by  using  die  seven  time  depeuient  variables  describing  the  motion,  and  the 
parameters  of  the  dual-yin  system.  These  additional  dynamical  quantities  are  stored  as  one 
dimensional  vectors  in  tile  copyuter. 

2.  Angular  Momentum 

Angular  momentum  was  derived  in  Equations  (123),  (124),  and  (125).  The 
angular  nxmaentum  is  a  vector,  but  for  verification  of  the  conservation  of  angular 
moinentum,  only  the  inagnitude  is  required.  Therefore,  the  three  vectorial  comptmems  for 
both  the  rotor  and  tiie  platform  are  conyuted  individuaUy  as  Hi,  H2,  H3,  Hip,  H2p,  H3p, 
and  then  the  magnitude  of  tiie  angular  momentum,  hQ|,  is  determined  and  stored. 

3.  Nutation  Angle 

The  nutation  angle  is  detennined  using  tiie  previously  determined  relation 


Since  it  is  a  function  ai  angular  momentum,  it  can  now  be  calculated  and  stored  as 

theta  []]. 

4.  Energy 

The  total  energy  of  the  system.  Equation  (136),  is  detennined  by  the  sum  of  the 
kinetic  energy,  as  derived  in  Equation  (134),  and  die  potential  energy  of  tiie  particle  masses 
of  die  mass-spiing-dashpot  danyer  system.  Equation  (135).  It  is  writtoi  as 

EMi-E  +  U-E+i/trJ  +  i*'/  (136) 

In  the  conyuter  program  the  kinetic  energy,  ke  [Q,  is  the  sum  of  die  constituents  of 
Equation  (134),  identified  as  T1  through  T13.  The  potential  energy  is  not  explicitly  stated 
in  the  program,  but  instead  is  included  in  the  calculation  ci  the  total  energy,  etotal  Q]. 
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5.  Core  Esergy  and  Poetulated  Core  &iergy 

The  platform  and  rotor  core  energy  of  die  dual-spin  system  is  defined  in 
Equations  (47)  and  (48)  respectively.  All  the  variaUes  are  availaUe,  such  dut  in  the 
computer  program  die  energies  are  readily  calculated  as  EcpQI  and  EcQl.  The  postulated 
nutadon  angle  may  now  be  computed.  Equation  (101)  indicates  a  sign  must  be  selected 
dqiending  on  whether  the  system  is  stable  or  unstaUe.  The  calculadon  of  nutation  angle 
over  time  in  Equation  (34)  will  indicate  stability.  The  computer  program  incorporates 
conditional  statements  to  assign  the  conect  sign  in  the  postulated  nutation  angle  eiqxession. 
EquaticMis  (101)  and  (106).  In  die  derivation  leading  up  to  the  postulated  nutation  angle,  o 
has  been  defined  as  the  relative  rotation  rate  of  the  rotor  with  reqiect  to  the  platform,  and  is 
a  positive  value.  This  provides  a  positive  contributim  to  the  angular  mcuiientum  vector, 
thus  providing  stability  to  the  systeoL  In  Mingoti’s  equations  of  motion,  the  reference 
coordinate  axes  are  fixed  to  die  rotor,  resulting  in  the  relative  rotation  of  the  platform  in  the 
counter-clockwise,  or  negative,  direction.  The  postulated  nutation  angle  equations  were 
developed  using  the  former  reference  firame.  To  compensate  for  the  difference  in  the 
reference  frames,  the  postulated  nutation  angle  equatitms  in  the  computer  program  have  o 
replaced  by -o. 

The  postulated  cote  energy  as  a  function  of  time  is  computed  for  the  platform  and 
the  rotor,  using  Equations  (140)  and  (141)  for  the  exponential  model  and  Equations  (142) 
and  (143)  for  the  Verhulst  model.  A  postulated  cme  energy  model  that  accurately  models 
dw  actual  core  energy  required  several  iterations  to  find  the  proper  values  d  the  exponential 
factor. 

Once  the  postulated  core  energy  models  are  computed,  the  Q'  and  Q  values,  as 
defined  in  Equations  (100)  and  (107),  are  determined.  The  postulated  nutation  angles. 
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ftah  [Q,  and  ata  (Q,  are  computed  using  Equation  ^101)  and  (106)  with  the  proper  sign 
previously  detetmined. 

The  revised  stability  criterion.  Equation  ^),  requires  the  time  rate  of  change  of 
the  core  energy  and  the  nutation  frequency.  Since  both  of  diese  parameters  may  vary  widi 
dme,  an  estimate  is  made  to  determine  if  die  modified  stability  criteria  conecdy  predicts 
stability  or  instability.  The  time  rate  of  change  of  the  core  energy  is  qiptmcimated  by  taking 
the  difference  of  the  final  and  the  initial  core  energy  values,  and  dividing  die  time  of  the 
simulation.  This  quantity  is  then  divided  by  dw  initial  nutation  frequency.  This  stability 
quantity  is  conqiuted  for  both  the  rotor  and  die  core,  but  is  indicated  in  Table  (1),  Summary 
of  Analyzed  Cases,  simply  as  a  negative  or  positive  value. 

B.  COMPUTER  PROGRAM 

Essential  to  the  verification  of  the  stability  criteria  is  the  computer  program  that 
integrates  the  dual-spin  system  equations  of  motion,  calculates  other  dynamical  quantities 
of  motion,  and  griqihs  the  results.  The  system  used  was  a  Sun  SPARC  2  workstation  with 
the  computer  code  written  in  C.  Intermediate  graphics  results  were  created  using  an  in- 
house  c(»nputer  graphics  program.  Final  graphics  output  was  performed  by  sending 
output  data  to  the  Deltagraph  graphics  package.  The  sequence  of  steps  of  the  computer 
code  are  explained  below.  The  con^iuter  code  is  included  as  Appendix  B. 

1.  Initialization 

The  main  computer  program  is  conpiled,  along  with  the  header  file  'rkk.h'  and  the 
function  'derivs*  immediately  preceding  it.  The  function  *derivs'  contains  Mingori’s 
equations  of  motion  rewritten  into  Equation  (145),  suitable  for  numerical  integration. 
Included  in  the  header  file  'rt(k.h'  is  the  numerical  integration  routine  'rk4,'  the  adrqrtive  step 
size  faction  'rkqc,'  and  the  driver  for  the  numerical  integration  routine  'odeint.'  Also 
included  are  functions  fen-  creating  and  freeing  the  vectors  and  matrices  used  by  the 
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oooqHiter  to  store  die  data,  and  a  functioo  for  enor  messages, 
e.  VuHabUs 

The  variables  required  for  both  the  numerical  integration  routine  and 
subseqiKntcalculatkms  are  defined  as  global  variables  before  the  main  program.  Variables 
used  only  within  a  specific  function  or  only  in  die  main  program  are  defined  in  their 
respective  funcdons.  Symbolic  constants  are  also  defined  before  the  main  program  using 
the  idefine  statement 

The  time  dependent  variables  of  motion  of  the  Mingori  system, 
o>i<  2,  i,  2\  i\  are  stored  in  a  seven  column  matrix.  All  other  parameters  to  be 

graphed  are  stored  in  one  dimensional  vectors.  Sunage  in  the  vectors  and  the  matrix 
permits  retrieval  by  the  graphics  subroutine  for  plotting  the  variable  over  time. 
b.  Input  FiU 

The  main  conqiuter  program  scans  die  input  file  fOT  the  necessary  parameters. 
AH  parameters  of  the  Mingori  system,  /i,  /i,  /s,  M,  m,  a,  *,  c,  l\\  /s',  M\m ',  a\  k\ 

c',  and  L,  are  controlled  with  the  input  file.  Also,  the  initial  conditions  for  the  time 
dependent  variables  of  the  system,  0)i,  a)2>  Ob,  Zt  z,  z',  z',  are  specified.  The  length  of 
time  of  the  simulation  and  a  variable  determining  the  desired  accuracy  are  specified. 
Finally,  the  exponential  factor  and  die  final  enogy  for  the  postulated  con  energy  functions 
are  read. 

2.  Preliminary  Calculations 

After  the  input  values  are  read  in,  initial  calculations  are  performed  prior  to  the 
numerical  integration  routine.  All  of  the  defined  terms  used  by  D.  L.  Mingori  [Ref  3]  in  his 
nmi-linear  equatitxis  of  motion  are  conqiured,  as  well  as  definititms  required  for  the  core 
energy  calculations,  /,,  //,  /,  total.  It,  W,  h  total,  and  hrigu. 

An  option  is  available  to  specify  critical  damping  in  either  the  platfmm  or  the 
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rotor.  By  using  the  number  1000.0  in  the  input  file  for  the  danqnng  coefficient,  the  main 
program  will  autonoatically  confute  the  coefficient  required  for  critical  danqiing  of  the 
mass-qmng-daslqxx  system  and  then  use  this  value  in  aU  subsequent  caknilarions. 

Also  computed  is  the  factor  'dxsav,  used  in  detennining  when  to  save  data.  The 
stq>s  between  evaluating  the  equations  of  moticm  can  become  small,  particularly  when  high 
accuracy  is  desired.  The  interval  required  for  gnq)hics  resolution  is  not  as  restrictive. 
Accordingly,  the  variables  are  »ved  only  if  the  step  is  greater  than  die  previously  saved 
step  by  the  factor  'dxsav.' 

3.  Numerical  Integration 

Mingoii's  nonlinear  equations  of  motion  are  numerically  integrated  by  the  fourth 
order  Runge-Kutta  method,  with  adaptive  step  size  cmitrol.  The  computer  code  used  is 
based  on  the  Runge-Kutta  method  listed  in  Numerical  Recipes  in  C,  [Ref  8].  The  adaptive 
step  size  control  permits  larger  integration  steps  during  smooth,  well  behaved  portions  of 
the  functions,  and  smaller  steps  during  the  more  irregular  sections  of  the  functions.  The 
integration  routine  accuracy  can  be  contrcdled  by  a  variable  in  the  input  file. 

Mingoii's  non-linear  equations  of  motion  are  contained  in  the  function  'derivs.' 
As  described  previously,  the  equations  have  been  rewritten  as  a  series  of  first  order  coupled 
differaitial  equations  suitable  for  numerical  integration. 

The  output  of  the  time  dependent  variables  of  die  system,  a>i.  o>z,  z,  i,  z\  i\ 
is  stored  in  an  array  of  seven  columns,  with  the  number  of  rows  required  a  function  of  the 
specified  time  interval  and  accuracy.  A  one  dimensional  vector  is  also  created,  with  the 
time  stored  for  each  step  saved.  This  permits  plotting  the  time  dependent  variables  and 
other  quantities  versus  the  time. 

4.  Calculation  of  System  Parameters 

Widi  the  array  of  the  time  depoidoit  variables  of  motion  as  a  function  of  time,  the 
other  system  quantities  listed  in  Section  A  may  now  be  calculated,  widi  the  values  stored  in 
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vecton  suitable  for  paphing. 

5.  Graphics  Output 

The  remainder  of  die  prognun  is  die  necessary  code  for  the  in-house  graphics 
program.  The  grqihics  program  is  initialized,  and  the  gnqiliics  window  is  opened.  The 
following  parameters  are  then  plotted  over  time:  n>i,  o>2,  n>3  -  m3  uiitiait  2,  i,  z',  i\ 

A  —  huudah  E  ~  Euutiali  ^ total  ~  ^total  iaitialt  initial*  pottalatad  ~~  pottalau^  initial* 

inMial  *  poatulata^  *  Ec initial*  ^*  V*  s»d  tf.  The  window  is  then  closed 
and  the  initial  conditions  and  peitinent  time  dqiendent  variables  and  other  quantities  of  die 
simulation  are  then  printed.  This  displays  the  gnqihical  representation  of  the  behaviOT  of 
the  dual-spin  system,  and  provides  the  actual  initial  and  final  values  of  the  variables  and 
associated  quantities. 

For  final  gr^ihics  output,  the  Deltagrtqih  gr^hics  program  is  utilized.  Conqiuter 
simulation  data  is  sent  as  an  ouqiut  file.  The  data  is  then  manipulated  into  a  suitaUe  graphic 
with  proper  scaling  and  axes  limits  to  best  represent  the  dynamics  of  the  computer 
siiiiulaticMi.  The  graphs  are  contained  in  Appendix  A. 

6.  Computer  Program  Validation 

The  two  aspects  of  the  cotrqiuter  code  requiring  validation  were  the  numerical 
integration  routine  and  MingOTi’s  equations  of  motion  with  its  associated  dynamical 
quantities  for  the  dual-qiin  system. 

The  validation  of  the  uamerical  integration  routine  was  perfOTmed  by  using  the 
equations  of  motion  for  a  simple  torque-fiee  axisymmetric  body.  Various  initial  conditions 
were  read  in,  and  the  dynamical  quantities  were  then  plotted  versus  time.  The  plotted 
behavior  was  then  con^iared  to  the  actual  response  of  the  system. 

The  validation  of  Mingoti’s  equadtms  of  motion  and  the  associated  dynamical 
quantities  required  a  mme  comprehensive  approach.  The  equations  of  motion  required 
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valfciatlon  for  proper  bduwMT  of  eaditimeHiependentviiu^  Then  die  kinetic 

enei:i^,  total  enetsy,  and  angular  momentum  must  be  verified.  Several  cases  were  run 
where  the  platfixm  mass  and  inertia  would  become  infinitesimally  smalL  Tlien  cases  were 
run  where  the  rotor  mass  and  inertia  would  become  infinitesimally  small  In  each  these 
cases,  subcases  were  run  where  the  mass-spring-dashpot  system  would  be  large, 
dominating  the  dynamics,  to  die  subcase  where  it  became  very  small,  so  die  system 
approximates  a  rigid  body.  These  different  scenarios  would  uncouple  and  isolate  the 
various  parts  of  the  equations  of  motion  to  verify  proper  derivation  of  die  equations.  In 
general,  the  platform  and  the  rotor  were  tested  under  ccmditions  varying  fimn  those  in  a 
rigid  body  scenario  to  those  in  a  lightly  damped  body.  The  system  was  dien  rested  as  a 
rigid  dual-^in  system,  as  a  dual-spin  syston  with  rotor  danqiing  only,  and  as  a  dual  spin 
system  with  platfonn  danqiing  only.  The  dual-spin  system  was  then  tested  with  both  the 
rotor  and  the  platfdxm  containing  dampers.  For  each  case,  the  time  dependent  variables  of 
motion  were  plotted  and  analyzed.  In  all  cases,  die  angular  momentum  was  compared  with 
the  initial  angular  momentum.  Since  all  cases  were  torque  free,  the  angular  momentum 
must  remain  constant.  In  all  cases  explored,  the  angular  momentum  verified  the 
correctness  of  the  equations  of  motion.  The  code  validation  cases  were  not  included  in  this 
thesis  due  to  the  large  number  of  graphs  and  data  that  were  required  to  establish  validation. 
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V.  ANALYSIS 


A.  INTRODUCTION 

1.  Objective 

The  devdopment  of  the  fevised  enosy-sink  stability  dwoiy  has  been  presented  in 
Chapter  n,  and  the  equations  of  motion  for  a  dual>qnn  system  are  contained  in  CSiapter  ni. 

die  numerical  integration  code  of  QuqMer  IV.  die  revised  energy  sink  staldity  themy 
can  now  be  verified.  The  core  energy  of  die  system  is  plotted  as  a  function  of  time  and 
conqnred  to  the  total  energy  for  agreement  with  theory.  The  staUlity  criterion  computed 
from  the  numerical  simuladon  must  then  agree  with  Equation  (97).  Conservation  of 
angular  momentum  is  used  in  each  case  to  verify  the  coirecmess  of  the  equations,  and  to 
ensure  the  accuracy  of  the  numerical  integraticMi  routine.  To  approximate  the  core  energies 
over  time,  an  exponendal  model  and  a  model  based  on  logistic  growth  are  then  explcved. 
The  logistic  growth  model  uses  an  equatimi  presented  by  Veihulst  [Ref  7],  and  is  referred 
to  as  the  Verhulst  model.  Postulated  models  of  core  energy  and  the  associated  nutation 
angles  are  dien  conqiared  widi  the  actual  core  oietgies  and  nutation  angle  for  agreement 

2.  Numerical  Simulation  Cases 

Four  distinct  cases  needed  to  be  addressed.  With  the  inertia  ratio  of  the  dual-spin 
system  greater  than  one,  a  stable  case  and  an  unstable  case  are  analyzed  (Cases  (1)  and 
(4)).  \ifith  the  inertia  ratio  less  than  one,  a  stable  case  and  an  unstable  case  are  also 
analyzed  (Cases  (2)  and  (3)).  For  these  four  cases,  an  exptmential  model  is  used  to 
postulate  the  core  energy  and  the  coireqxniding  nutation  angle.  For  the  case  of  the  inertia 
ratio  greater  than  one  and  unstable,  and  the  two  cases  of  the  inertia  ratio  less  than  one, 
stable  and  unstable,  the  Verhulst  iiKxlel  is  postulated  (Cases  (S),  (6),  and  (7)).  The  case  of 
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die  inertia  ntio  greater  dian  cne  and  stabte  was  not  included;  from  the  otceUem  agreement 
of  Cases  (5).  (6).  and  (7),  one  can  see  diat  excellent  agreement  would  also  occur  for  this 
trivial  case,  making  it  unnecessary  10  include.  All  seven  cases  oe  summariaed  in  TaUe  (1). 
For  each  case,  die  dynamical  quantities  are  plotted  versus  dine  until  the  system  reaches  a 
stable  state.  Each  quantity  is  then  plotted  for  the  first  one  hundred  seconds  to  show  detail. 


CASE 

h  total 
htotal 

STABILITY 

X 

X' 

AiODEL 

COMMENTS 

1 

1.039 

negative 

negative 

exponential 

large  danqnng 
inplatfonn 

2 

0.907 

stable 

n^ative 

negative 

exponential 

large  damping 
in  platform 

3 

0.897 

unstable 

positive 

positive 

exponential 

insufilcient  dancing 
in  platform 

4 

1.025 

unstable 

positive 

podtive 

esqxmential 

insufficient  dandling 
in  platform 

5 

0.907 

staUe 

negative 

negative 

Veifaulst 

same  conditions 
asCase2 

6 

0.897 

unstable 

positive 

positive 

Verhulst 

same  conditions 
asCase3 

7 

1.025 

unstable 

positive 

positive 

Verhulst 

same  conditions 
asCase4 

Table  (1)  Summary  of  Analyzed  Cases 


In  Appendix  A,  (2)  through  (8)  list  the  system  parameters,  core  energy 
parameters,  and  initial  condiuons.  immediately  flawing  these  tables  are  the  gnyihs  of  the 
inqiortant  dynamical  quantities.  A  typical  geosynchronous  dual-spin  satellite  was  selected 
for  the  numerical  simulation.  The  platfcxm  and  rotor  masses  and  inertias  are  listed  as  part 
of  the  initial  conditicms,  and  are  the  same  for  all  cases.  The  inertia  ratio  is  made  greater 
than  or  less  than  one  by  selecting  L,  die  distance  between  die  rotor  and  platform  centers  of 
mass,  to  be  0.3  or  1.0  meter  respectively.  As  a  default  set  of  values,  all  of  the  mass- 


qiniig«dBdipot  synems  htvt  m  and  m*  equal  to  1.0  kg.  A  and  k'  equal  to  1.0  NAn,  and  c 
and  coequal  to  1.0  kgfwec.  To  establish  the  stable  cases,  additional  energy  dissipation  is 
required  in  the  platform.  To  achieve  diis.  cases  (1),  (2).  and  (5)  have  m'  increased  to  20.0 
kg  and  tf  increased  to  10.0  kg/sec.  In  all  seven  cases,  die  initial  conditions  are  die  same. 
The  platform  rotates  at  a  geosynchronous  angular  velocity,  the  rotor  qnns  at  the  higher  rate 
ci  1.5  rad/sec.  and  a  perturbation  is  introduced  by  an  initial  transverse  angular  velocity  of 
0.1  rad/sec.  The  mass>spring-dashpot  systems  have  no  inidal  displacement  or  initial 
velocity.  The  core  energy  parameters  are  also  listed  in  the  taUes.  The  initial  core  energies 
are  determined  tty  the  system’s  initial  conditions.  The  final  core  energies  were  taken  from 
the  numerical  simulation  data.  The  exponential  factors,  r  and  were  then  determined 
through  an  iterative  process  to  best  fit  the  modeled  core  energy  to  die  actual  core  energy. 

B.  DISCUSSION 

The  individual  cases  can  now  be  analyzed.  The  graphs  of  the  dynamical  quantities 
are  explored,  and  the  data  will  affirm  die  revised  stability  thecny. 

1.  Angular  Momentum 

For  each  case,  the  angular  momentum  is  plotted  versus  time.  Figures  (5),  (6), 
(16).  (17).  (27),  (28),  (38),  (39),  (49),  (50),  (60),  (61),  aD.  and  (72).  Because  the  dual- 
^in  system  has  no  external  forces,  angular  momentum  must  be  ccxiserved.  For  the  stable 
cases,  there  is  excellent  agreement,  with  angular  momentum  varying  by  less  than  one  one- 
hundredth  of  a  percent  over  the  length  of  the  data  run.  This  confirms  the  equations  of 
motion  and  validates  the  accuracy  of  the  numerical  integration  routine.  For  the  unstable 
cases,  the  angular  momentum  percent  difference  increases  to  qiproximately  four  one- 
hundredths  of  a  percent  This  is  attributed  tt>  the  increased  dynamics  of  the  system  as  it 
establishes  the  qiin  about  die  transverse  axis,  introducing  veiy  small  errors  in  the  numerical 
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imegmioii  routiiie.  The  enon  remain  voy  mall,  and  the  equatioas  of  motion  and  die 

iMiiwiitirMl  itiiiyAni  wiwrfti#  wwwtw  anw—li* 

2.  Total  Energy,  Platform  and  Rotor  Core  Energy 

The  gnqihs  eaagy  are  in  Hgures  (7)  throat  (10).  (18)  through  (21).  (29) 
through  (32).  (40)  through  (43).  (SI)  through  (54).  (62)  through  (65).  and  (73)  through 
(76).  Several  obsorvations  can  be  made. 

The  total  energy  curve  reaches  equilifarium  before  bodi  die  {datfonn  and  die  rotor 
core  energies  do.  For  Case  2  and  5.  die  total  energy  iqipears  to  reach  a  maximum  and  then 
drop  down  to  a  final  equilibrium  value.  Careful  conqiaxiscMi  of  Figures  (18)  and  (49), 
shows  the  same  system  with  the  same  initial  conditions,  but  with  a  slighdy  differait  curve. 
It  can  be  deduced  diat  the  energy  is  not  steady,  but  is  still  oscillating.  The  sampling 
frequency  coincidentally  saved  values  near  the  same  magnitude  ci  energy  in  the  region 
from  about  3000  seconds  to  10000  seconds. 

The  total  energy  for  C4ses  (1).  (3).  (5).  and  (^  decreased.  For  Cases  (2),  (4), 
and  (7),  rqaesenting  both  staUe  and  unstable  systons,  it  increased.  This  can  be  explained 
easiest  with  the  use  of  equations  (26)  and  (27).  For  Crises  (2),  (4),  and  (7),  the  system  is 
settling  out  about  the  axis  with  the  minimum  moment  oi  inertia.  Because  there  are  no 
external  forces,  the  angular  momentum  is  constant  In  die  equadon  for  angular  momentum. 
Equation  (26),  the  inertia  terms  are  squared.  But  in  the  equation  for  energy.  Equation  (27), 
the  inertia  terms  are  not  squared.  Therefore,  as  die  angular  velocity  transfers  to  the  axis 
with  the  tninimiiin  moment  of  inertia.  Equations  (26)  and  (27)  show  that  die  energy  will 
increase.  These  equations  apply  to  a  sinqile  dual-spin  system.  Equations  (123)  and  (136) 
for  the  Mingori  dual-spin  system  would  show  the  same  result,  provided  the  energy 
absorbed  by  the  mass-qiring-dashpot  system  and  the  energy  associated  with  the  motor  is 
tess  dian  the  energy  increase  associated  with  transferring  die  spin  to  the  axis  of  minimum 
moment  of  inertia.  This  is  the  situation  for  Cases  (2),  (4),  and  (7). 
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The  total  energy  value  will  always  be  between  diat  of  die  platform  core  energy 
and  the  rotor  cme  energy.  The  platform  core  energy  will  have  smaller  values  as  its 
equation  assumes  dut  the  rotor  is  rotating  at  the  same  rate  as  the  platform,  diereby  not 
accoundngforthelargeamountofkineticenergy  associated  with  the  rotor.  The  rotm  core 
energy  will  be  higher  than  the  total  energy  as  it  assumes  that  the  slowly  spinning  platfmm 
is  rotating  at  the  same  rate  as  the  rotor,  providing  die  system  with  additional  kinetic  energy. 

The  hundred-seconds  graphs  d  total  energy  versus  time  shows  curves  with 
several  different  behaviors.  The  time  rate  d  change  of  total  energy  includes  the  energy 
dissipation  rate  of  the  rotor  and  platform  mass-spring-dashpot  systems,  and  the  rate  of 
work  due  to  the  motor  torque  maintaining  the  relative  rotadon  rate.  In  the  hundred-seconds 
gnqihs,  aside  frcxn  the  general  slc^  of  the  curve,  diere  is  no  apparent  coirelation  between 
the  specific  behavior  of  the  total  energy  to  that  of  the  core  energies.  Any  relationship  that 
exists  is  masked  by  the  motcu’  and  mass-spring-dashpot  system’s  influences  on  total 
energy. 

3.  Stability  Criterion 

The  revised  stability  criterion  of  equation  (97)  states  that  die  dme  rate  d  change  of 
the  core  energy  over  the  respective  nutadon  frequency  must  be  less  than  or  equal  to  zero. 
The  sign  of  die  stability  criterion  fOT  each  case  must  be  determined.  The  dme  rate  of  change 
of  the  core  energy  was  determined  by  dividing  the  final  less  the  inidal  value  of  the  core 
energy  by  the  length  of  dme  of  the  case.  The  nutadon  frequencies  were  computed  using 
initial  conditions.  The  sign  of  the  stability  criterion  is  listed  in  Tables  (1)  through  (8).  The 
numerical  value  was  not  listed  because  the  above  assumptitxis  used  during  the  computation 
give  it  no  merit  Cases  with  inertia  ratios  very  close  to  one  were  intentionally  selected  to 
test  the  inertia  ratio  in  this  transition  region,  bi  all  cases  analyzed,  the  sign  of  the  revised 
stability  criterion  was  consistent  widi  die  stability  of  the  systeia 
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4.  Postulated  Cmv  Energy  and  Nutation  Angle 

An  exfdidt  relationship  for  core  enagy  as  a  functioa  of  tune  does  not  exist  If  an 
equation  describing  core  energy  did  exist,  then  by  using  Equations  (101)  and  (106)  the 
nutation  angle  as  a  function  of  core  energy  and  time  could  be  predicted  for  a  dual-spin 
system.  This  could  then  be  extended  to  an  actual  dual-spin  satellite.  Given  a  sufficient 
modd  die  satellite,  the  staldity  and  the  nutadon  angle  as  a  function  of  time  could  be 
predicted.  An  objective  of  diis  thesis  is  to  see  if  an  equation  for  the  core  energy  can  be 
develqied  that  would  adequately  describe  die  nutation  angle  as  a  functitMi  of  time  for  both 
the  stable  and  unstable  condidons.  The  eiqxmendal  model  and  the  Verhulst  Iqgisdc  model 
were  explored. 

4.  Expotuntial  Core  Energy  Model 

Observadon  of  the  core  energy  as  a  funetkm  of  time  for  a  stable  system  would 
lead  one  to  ctxiclude  that  it  bduives  in  an  exponential  manner.  Equadons  (140)  and  (141) 
are  exponendal  representadems  of  rotor  and  platform  energy  as  a  funedon  of  dme.  The 
initial  core  energy  was  determined  by  the  initial  condidons  of  the  dual-spin  system.  The 
final  core  energy  was  determined  by  tunning  die  numerical  simuladon  and  calculating  it  at 
the  end  of  the  simuladon.  If  this  was  not  available,  a  value  could  be  estimated.  Applying 
the  principle  of  conservadon  of  angular  momentum  and  noting  that  the  system  will 
eventually  spin  about  one  of  the  primary  axes,  then  the  equadtms  for  angular  momentum 
and  total  energy  can  be  used  to  solve  for  the  angular  velocity  about  that  axis.  Subsdtudng 
this  into  Equation  (47)  <x  (48),  one  would  arrive  at  an  estimated  final  core  energy. 
Although  it  would  not  be  die  actual  final  core  energy,  as  motor  torque  contribution  and  the 
mass-spring-dashpot  system’s  energy  dissipation  was  not  accounted  for,  it  would  be 
sufficiendy  close  to  satisfy  the  computational  requirements.  The  exponential  factor  is 
dependent  upon  the  system  parameters  and  initial  conditions.  For  the  cases  presented, 
different  values  were  tried  until  good  agreement  was  establiriied  with  the  core  energy 
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curve.  Additkaalieieiich  woukibeiequiicd»detenniiicasuiiabieejq>oncotialfaciorfor 
an  actual  laldlite,  taking  into  account  die  parameters  and  die  initial  conditions.  Cases  1 
dnou^  4  contain  die  numerical  siinulatian  data  for  the  eiqxxiential  model  Theactualand 
postulated  curves  for  core  energy  are  contained  in  Figures  (1 1)  through  (IS),  (22)  duough 
(24),  (33)  through  (37),  and  (44)  through  (48).  For  cases  1  and  2,  both  stable,  diere  is 
excellem  agreement  between  die  core  energy  and  the  postulated  core  energy.  This  in  turn, 
results  in  excellent  agreement  between  the  actual  nutation  angle  and  the  modeled  nutation 
angle.  The  actual  nutation  angle  is  determined  using  the  angular  momentum  quantities,  as 
shown  .  Equation  (34).  The  time  dependent  variable  required  to  cooqiute  the  nutation 
angle  is  the  angular  velocity  about  the  spin  axis.  The  modeled  nutation  angle  is  determined 
in  Equations  (101)  and  (106),  and  is  a  fiinctimi  of  only  one  time-dependent  variable,  core 
energy.  Herein  lies  the  potential  of  the  core  energy  theory.  A  sufficient  model  of  core 
energy  over  time,  as  in  Cases  1  and  2,  will  provide  an  excellent  prediction  of  nutation 
angle,  without  requiring  any  knowledge  the  specific  angular  velocities  of  the  system  as  a 
function  of  time. 

Cases  3  and  4  illustrate  the  exponential  model  for  an  unstable  dual-spin 
system.  The  exponential  model  for  cott  energy,  and  its  associated  nutation  angle,  rapidly 
approach  their  final  values.  The  actual  core  energy  and  nutation  angle,  however,  behave 
quite  differently.  The  initial  conditions  have  the  system  near  an  unstable  equililxrium.  The 
system  moves  from  the  unstable  to  the  stable  equilibrium  slowly  at  first.  It  then  increases 
the  rate  at  which  it  iq>proaches  stable  equilibrium,  passes  through  an  inflection  point,  and 
then  approaches  equilibrium  asymptotically.  It  is  clear  that  the  exponential  model 
represents  this  behavior  poorly.  Verhulst’s  logistic  equation  was  then  addressed  to 
determine  its  adequacy  in  modeling  the  dual  spin  system. 
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k.  Vtrkmbt  iMgisde  C»n  Energy  Model 


The  Veilittbtk^isdceqiiiation  was  introduced  in  ChqMer  in.  Hgure(4)and 
die  assodaieddesciqxion  explains  the  diaracieristics  of  die  equation.  The  genoal  sluqie 
the  curves  in  Figure  (4)  is  very  similar  to  the  staUe  and  unstaUe  cases  of  the  dual>^in 
system.  Case  S  is  die  Veihulst  model  of  die  Case  2  stable  system.  Hguies  (SS)  throu^ 
(59)  show  diat  the  Verhulst  logisdc  equation  can  achieve  ocellem  agreement  widi  core 
energy  and  nutadon  angle,  just  as  die  exponential  model  did. 

Cases  6  and  7  are  the  same  as  Cases  3  and  4,  except  the  Verhulst  logisdc 
equadon  is  used  to  model  the  core  energies.  Figures  (66)  dirou^  (70)  and  (77)  through 
(81)  illustrate  die  modeled  and  actual  core  energies.  Although  the  Verhulst  model  for  rcxor 
core  energy  had  excellent  agreement,  it  performs  poorly  when  modeling  the  unstable 
system.  The  rotor  core  energy  begins  at  an  initial  value  and  then  decreases  to  an 
equilibrium  value.  Referencing  Figure  (4),  the  Verhulst  logisdc  equadon  will  model  the 
rotor  core  energy  exponentially. 

To  take  advantage  of  die  Verhulst  logisdc  equadon’s  curve  beginning  near  the 
the  unstable  equililxium  and  its  progressicm  to  the  asynqitodcally  stable  equilibrium,  the 
con  energy  must  increase  over  dme.  The  platfcmn  core  energy  behaves  in  this  manner  as 
the  nutadon  angle  goes  firom  a  small  angle  to  ninety  degrees.  Figures  (68)  through  (70) 
and  (79)  through  (81)  show  the  Verhulst  modeled  and  the  actual  platform  core  energies, 
and  the  modeled  and  the  actual  nutadon  angles.  The  agreement  between  the  actual  and 
modeled  energies  and  nutadon  angles  was  very  good.  The  general  sh^  of  curve  was 
ccmsistent,  with  only  a  slight  deviadcm  in  the  center  of  the  curve,  and  then  a  small  deviadon 
as  tte  nutation  angle  approaches  the  equilibrium  value  of  ninety  degrees.  With  this 
agreement,  it  is  established  that  a  core  energy  model  exists  that  can  represent  both  stable 
and  unstable  cases.  For  each  of  the  cases,  con  energy  is  plotted  versus  nutadon  die  angle. 
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Rgufes  (15),  (26),  (37),  (48),  (59),  (70),  and  (81).  The  gnphici  package  pennitted  only  a 
semi-log  plot  of  rotor  core  eneigy.  The  plot  did  not  provide  any  insight  on  rdadonships 
for  die  dynamics  of  the  dual-spin  system.  For  both  the  stable  and  unstable  cases,  die  log- 
log  pltH  between  platform  core  eneigy  and  the  nutation  angle  is  linear  from  ninety  degrees 
until  about  two  degrees.  For  nutation  angles  less  than  two  degrees,  the  platform  core 
energy  asynqitotically  approaches  the  equilibrium  value.  Case  5  and  6  is  the  same  system 
widi  the  same  initial  conditions,  with  the  excepticm  of  die  increased  damping  system  in  the 
platfmm  for  C^tse  5.  The  initial  platform  core  energies  and  nutation  angles  are  very  nearly 
the  same  for  each  case,  as  shown  in  Figures  (57),  (58),  (59),  (68),  (69),  and  (70).  One 
can  conclude  by  this  observation  and  the  definitkm  of  nutation  angle  as  a  function  of  core 
energy.  Equation  (101),  that  there  exists  a  continuous  platform  core  energy  versus  nutation 
angle  curve.  The  appearance  would  look  similar  to  the  curve  created  by  splicing  Figures 
(59)  and  (70)  together.  By  varying  CMie  parameter,  for  example  platform  eneigy  dissipation 
rate,  the  system  would  progress  along  this  curve  and  achieve  equilibrium  with  a  zero 
degree  nutation  angle  or  achieve  equilibrium  with  a  ninety  degree  nutation  angle.  This  is 
what  was  done  with  cases  5  and  6.  Adding  as  the  diird  dimension  to  die  curve  the  time  rate 
of  change  of  core  energy  would  then  reveal  the  stalnlity,  die  initial  direction,  and  the  rate  at 
which  the  system  will  arrive  at  the  equilibrium  condition. 

C.  FURTHER  RESEARCH 

The  revised  stability  criterion  for  a  dual-spin,  quasi-rigid,  axisymmetiic  system  was 
established.  Numerical  simulation  was  then  used  to  verify  the  revised  stability  theory. 
Further  research  could  be  conducted  in  several  areas.  The  specific  contributions  to  the  total 
eneigy  could  provide  some  insight  How  much  energy  and  in  what  manner  does  the  motor 
torque  contribute  to  the  total  energy  for  both  the  stable  and  the  unstable  cases?  Also,  by 
plotting  the  energy  dissipation  system’s  contributions  over  time,  one  could  determine  its 
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effect  on  total  energy,  as  well  as  conqtaring  it  to  the  cote  energy  theocy  in  Equation  (97). 

Relati(»ish4>s  for  cote  energy  as  a  fiincticm  of  dme  were  explored.  It  was  established 
that  the  Veihulst  logistic  equation  could  be  aj^lied  to  a  stable  dual-spin  systan.  By 
changing  its  parameters,  this  equaticHi  also  applied  to  the  unstable  dual-spin  system  with 
very  good  agreement  Further  research  could  be  directed  to  find  one  equation  that  would 
address  both  the  stable  and  unstable  cases  of  the  dual-spin  system  without  changing  its 
parameters.  The  equation  for  logistic  growth  with  a  direshold  [Ref.  7]  as  ^iplied  to  the 
platform  cme  energy  shows  promise  in  this  regard. 

ic  =  -r'il - - (146) 

\  Eciuutabulx  Ec  final  f 

Given  the  i  initial  conditions,  the  platform  core  energy  will  typically  start  at  some 
intermediate  vaiue  and  will  either  increase  or  decrease  as  the  dual-spin  system  reaches 
equilibrium  with  a  nutation  angle  of  either  zoo  degrees  or  ninety  degrees.  Equation  (146) 
would  be  well  suited  to  model  the  platform  core  energy.  There  exists  a  platform  core 
energy  value,  Eq  Hnsuibu\  where  the  system  is  at  an  unstable  equilibrium.  Equation  (146) 
shows  that  at  exactly  this  value,  the  rate  of  change  of  the  platform  core  energy  will  be  zero. 
For  any  value  below  the  unstable  equilibrium  value,  the  core  energy  will  approach  the  value 
of  zero.  For  Case  5,  the  final  core  energy  for  the  stable  condition  was  0.122  J.  Although 
this  final  condition  cannot  be  represented  in  Equation  (146),  it  is  sufficiently  close  that  the 
equation  may  still  be  used  to  represent  the  platfcmn  core  energy.  For  any  value  above  the 
unstable  equilibrium  value,  the  core  energy  wiU  approach  die  equilibrium  value  of  Ec  final* 
associated  with  the  nutation  angle  at  ninety  degrees.  Initial  conditions  would  determine  the 
initial  core  energy.  The  final  core  t  Tgy  can  be  estimated  as  described  earlier  in  this 
chapter.  Finally,  the  exponential  factor  /  may  then  be  determined  based  on  the  system 
parameters  and  initial  conditions. 

AddititMial  analysis  can  be  performed  on  the  revised  stability  criterion.  The  time  rate 
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of  duuife  of  core  energy  over  nutational  fiequency  may  be  calculated  and  plotted  versus 
time  or  versus  other  parameters.  The  behavior  of  diisstabili^  criterion  and  die  magnitude 
ai  it  for  various  condititxis  could  provide  some  insight  on  postulating  a  relationship  for 
core  energy  as  a  function  of  time. 
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VL  CONCLUSION 


The  exisdiig  enosy-sink  stability  criierion  was  inirodiiced  as 

(41) 

It  was  shown  that  an  inconsistaKy  in  the  develtqmient  disproves  the  existing 
assumption  Aat  die  motor  energy  input  exacdy  balances  shaft  frictional  losses.  A  revised 
energy-sink  stability  ciiterimi  was  then  developed  based  on  Hubert’s  definition  of  core 
energy  and  was  presented  as 

(97) 

This  criterion  compliments  die  existing  dieory.  Numerical  simuladon  was  required  to 
validate  the  theory.  The  Mingori  dual-^.  quasi-rigid.  axisymmetric  system  was  selected 
for  the  numerical  simuladon.  Several  cases  were  analyzed  to  verify  the  revised  energy-sink 
stability  criterion.  By  correcdy  postulating  the  platform  or  the  rotor  core  energy,  the 
stability  of  the  system  could  be  determined.  Specific  knowledge  of  the  energy  dissipadon 
rates  for  both  the  platform  and  the  rotor  are  no  longer  required. 

An  exponendal  model  and  the  Veihulst  logistic  model  for  core  energy,  and  their 
reladonship  to  the  nutadon  angle,  were  explored.  The  exptmendal  model  had  excellent 
agreement  with  the  staUe  cases,  but  was  inadeq******-  in  representing  die  unstable  cases.  The 
Veihulst  logistic  modd  established  diat  an  explicit  reladonship  for  core  energy  could  be 
devdt^ied.  Nutadon  angle  as  a  frincdon  of  core  energy  for  the  dual-^nn  system  could  then 
bepttdicied.  The  excellent  agreement  of  die  postulated  core  energy  and  nutadon  angle  with 
the  actual  core  energy  and  nutation  angle  confirms  the  revised  energy-sink  stability 
criterion.  Additional  research  is  required  to  find  an  optimum  equation  for  erne  energy  as  a 
function  of  time  to  rq;nesent  bodi  the  stable  and  unstable  cases. 
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APPENDDCA  NUMERICAL  SIMULATION  DATA 


System  Parameters 


Initial  Conditions 


Flatfonn 


Rotor 


Platform 


Rotor 


1600  kgm2 

/i- 

1000  kgm^ 

V*  1500  kgm2 

/3* 

1200  kgm2 

AT-  1000.0  kg 

700.0  kg 

m's  20.0  kg 

m  » 

1.0  kg 

a's  1.0  m 

aa 

1.0  m 

k'*  1.0  S 

m 

k» 

1.0  £ 
m 

C-  10.0^ 

c* 

1.0 

sec 

Ls  0.3  m 

total 
ff  foca/ 

=  1.039 

z'  « 

0.0  m 

z  » 

0.0  m 

z'  * 

0.0  f 

z  » 

0.0  m 

7.27x10-5^ 

1.5^ 

©1  * 

0.10  od- 

sec 

0.0rg 

Core  Energy  Parameters 


Platform 


Rotor 


£c  iiutial  ~ 

13.402  J 

initial  “ 

3145.4  J 

final  ~ 

0.072  J 

final  ® 

3161.7  J 

iSc"  _ 

A' 

negative 

II 

negative 

/  * 

-.00134  s-» 

r  = 

-.00134  s-> 

Table  (2)  Case  1  Parameters 


energy  (J)  I  |  energy  (J) 


Figure  (5)  Total  Enogy  and  Pooent  Difference  Angular  Momentum  Versus  Time 


Case  1:  Is/It>l,  Stable,  Exponential  Model 
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Figure(6)  Total  Eneigy  and  Percent  Difference  Angular  Momentum -Hrst  100  Seconds 
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angular  momentum  difference  (percent) _ 1  I _ angular  momentum  difference  ftieicent) 


«*»iy(J)  I  I  energy  (J) 
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Figure  (7)  Total  Enogy  and  Rotor  Core  Energy  Versus  Time 


Case  1:  Is/It>l,  Stable,  Exponential  Model 


—  Maleneigy 


•~-iolor  core  energy 


31483 


31473 


31463 


314S3 


Figure(8)  Total  Energy  and  Rotor  Core  Eneigy  First  100  Seconds 
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core  energy  (J) _ I  I _ core  eneigy  (J) 


enaiy(J)  I  metgyiJ) 


coieeiieffy(J) 


cofeeiiergy(J)  ]  I  cancneigy(J) 


Figure(ll)  Moddedand  Actual  Rotor  Core  Energy  and  Nutation  Angle  Versus  Tinx 


Figure  ( 12)  htodeled  and  Actual  Rotor  Core  Energy  and  Nutation  Angle*  First  100  Seconds 


72 


nutetko  angle  (degrees) _ I  I  iiul*ik]n  angle  (degraea) 


coie  energy  (J) 


Figure  ( 13)  Modeled  and  Actual  Platform  Core  Energy  and  Nutation  Angle  Versus  Time 
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Dutetkm  angle  (degieea) 


conoHffyCJ) 


coicenefsy(J) 


CASE  2:  ^  <  1  ,  Stable,  Exponential  Model 


System  Parameters 


Platform 


Rotor 


Initial  Coodiucos 


Platform 


Rotor 


/i'>  1600  kgm2 

/i« 

1000kgiii2 

V*  1500  kgm2 

/3* 

1200  kgm2 

1000.0  kg 

Afs 

700.0  kg 

20.0  kg 

m* 

1.0  kg 

o'  s  1.0  m 

as 

1.0  m 

*'=  ‘OS 

ks 

1.0  ^ 
in 

c'.  10.0^ 

c  = 

1.0 -Is*. 

sec 

1.0  m 

^atoui  ^ 

hwud 

0.907 

/  * 

0.0  m 

z  « 

0.0  m 

i'  » 

0.0  m 
s 

z  * 

0.0  m 

©s'* 

7.27  X  10^  nj 
see 

OJJ  * 

isod 

‘•^sec 

a>i  » 

O'OSJ 

0)2* 

0.0  od. 

sec 

Gore  Eneiigy  Parameters 


Platform 


Rotor 


Ec  initial  * 

15.34  J 

Ec  initial^ 

3147.3  J 

Ee/mta  = 

0.123  J 

Ecfinat  * 

3170.9  J 

negative 

£c 

negative 

■■■■■  S 

A' 

X 

r'  = 

-.00102  s-i 

r  s 

-.00102  s-» 

Case2ParB 

imeters 

Figure  (17)  Total  Energy  aiKl  Pocent  Difference  Angular  Momentum  •  Hrst  1()0  SecoiKls 
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energy  (J)  f  energy  (J) 


coie  energy  (J) 


it)SMan»  I  (r)XSjai» 


cofeenerfyCJ)  I  I  careeiwf|y  (J) 


con  energy  (J)  I  I  coraeaefgy(J) 


Fisure(22)  ModdedandActualRotorCoKEnenyandNutation  Ancle  Venus TiiiM 


31S0 

3149J 

3149 

31483 

3148 

31473 

3147 

31463 

3146 


Figure  (23)  Modded,  Actual  Rotor  Core  Energy  and  Nutation  Angle -First  100  Seconds 
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nutation  angle  (degrees)  I  ”  I  nutation  angle  (degrees) 


caieeiiHgy(J)  I  17)  I  concoeigy(J) 


nutation  ingle  (degrees)  _ I  n  I _  nutation  angle  (degrees) 


cape«nHfy(J) 


caieenHfy(^ 


'4 


CASES:  ^<1  ,Uii8tabit,Eipoiicntial Model 


Sysmn  PusinctErs 


Initial  Conditions 


Haifonn 


Rotor 


I\  »  1600  kgn^ 

h  * 

lOOOkgixP 

V-  lS00kgm2 

/3- 

1200  kgnP 

1000.0  kg 

700.0  kg 

1.0  kg 

m* 

1.0  kg 

o'*  1.0  m 

o* 

1.0  m 

H-  1.0  H 

10  S 

m 

c'*  LO  ^ 
sec 

c* 

1.0.^ 

sec 

1.0  m 

IzuuLs 

hiout 

0.897 

Platform 


Rotor 


Cote  Energy  Panuneteis 


Platform 


Rotor 


^emitial 

15.1  J 

3061.6  J 

Eefima  * 

1 160.9  J 

-  1489.9  J 

positive 

£c 

“  »  positive 

A 

K  * 

-.0005  s-» 

r  «  -.0005  S-* 

— 

Table  (4)  Case  3  Parameters 


eaeify(J) 


ngure(27)  Total  Enogy  and  PacentDifSaenoe  Angular  Momentum  Versus  Tim 


Flgure(28)  Total  Eneij^  and  Percent  Difierence  Angular  Momentum -Rrst  100  Seconds 
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angukrnioawnlitmdiflbfence(peic«l)  I  ”  I  aofubriiioiiiealiimiliflcnaoe(|MioeDt) 


Figure  (31)  Total  Energy  and  Platfmn  Core  Energy  Versus  Time 
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Rguie(32)  Total  Energy  and  Platform  Core  Energy  -  First  100  Seconds 
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Case  3:  Is/b<l,  Unstable,  Exponential  Model 
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igure(3S)  Modeled  and  Actual  Platform  Core  Energy  and  Nutation  Angle  Versus  Time 


Hgure  (36)  Modded,  Actual  Platfcnm  Core  Energy  and  Nutation  Angle-  First  1(X)  Seconds 
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nutaiton  angle  (degrees) 


caneae^nr(J) 


coieetiefgy(J) 


CASE  4:  ^  >  1  ,  Unstable,  Exponential  Model 


Table  (S)  Case  4  Parameters 


initial  Conditions 
Platfonn  Rotor 


i  * 

0.0  m 

z  * 

0.0  m 

i'  * 

00? 

i  * 

0.0  a 

S 

tori* 

7.27  X  10-*  gi 

0)j- 

‘“’sec 

CDi  B 

0.10 

sec 

a>a* 

0.0 

sec 

Gore  EneiQ' Parameters 
Platform  Rotor 


mj 

EciMal*  3059.7  J 

Ecfinal  - 

1247.1  J 

^c/ina/  1552.5  J 

positive 

~  »  positive 

X 

/  * 

-.0005  s-i 

r  »  -.0005  S"* 

System  Parameters 
Platform  Rotor 


/i'«  1600  kgm2 

/j  s 

1000  kgnP 

/s'*  1500  kgn^ 

/s  * 

1200  kgm2 

AT*  1000.0  kg 

M« 

700.0  kg 

m'»  1.0  kg 

m* 

1.0  kg 

o'  *  1.0  m 

a- 

1.0  m 

log 

10  S 

m 

c'.  1.0^ 
sec 

c* 

1.0  Jst 

sec 

L=  0.3  m 

fsteul  — 
h$0tal 

1.025 
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(f)XSjai»  _ I  I _  (r)XfjMia 


■ngultf  momentum  difference  (perocnt) _ I  I _ angular  momentum  difference  (perecol) 


energy  (J) 


Figure(40)  Total  Energy  and  Rotor  Core  Eneigy  Versus  Time 
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Figure  (41)  Total  Energy  and  Rotor  Core  Energy  -  F%st  100  Seotmds 
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Figure  (44)  Modeled  and  Actual  Rotor  Gxe  Energy  and  Nutation  Angle  Versus  Tiiro 
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Figure  (45)  Modded,  Actual  Rotor  Core  Energy  and  Nutation  Angle  •  First  1(X)  Seconds 
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coi«eiieify(J) 


System  Bminetm 
Plitfionn  Rotor 


v-  1600  kfm2 

/i- 

1000  kgnP 

V-  1500  kgm2 

/3- 

1200  kgo^ 

1000.0  kg 

700.0  kg 

ir'»  mo  kg 

m» 

1.0  kg 

1.0  m 

fl* 

1.0  m 

a*  1.0^ 

m 

k« 

10  S 

m 

c'*  10.0  ^ 
sec 

c« 

1-0^ 

sec 

1.0  m 

ftioiai 

Ifiatal 

«  0.907 

Tnirial 

Platform  Rotor 


/  - 

0.0  m 

z  » 

0.0  m 

r  - 

0.0  a 

0 

z  « 

0.0  ni 

S 

m3'* 

7.27  X  1(H  ad 

SOC 

mj* 

1.5  gd 

sec 

0.10  aj 

IK 

0.02^ 

sec 

Gore  Enasy  Panmetas 
Platform  Rotor 


^emUal  * 

15.341  J 

Ecimtiaf  3147.3  J 

EcfiMl  * 

0.122  J 

Ecyw  -  3170.9  J 

A' 

n^ative 

—  *  negative 

A 

/  « 

-.0002  s-» 

r  *  -.0002  s-» 

TaUe(6)  Case  S  Patameteis 
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!^gufe(49)  Total  Enogy  and  tacentDiffiereooe  Angular  Momentum  Venais  Time 
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Figure(50)  Total  Energy  and  Percent  Differenoe  Angular  Momentum -Rrst  100  Seconds 
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Figure(53)  Total  Eneigy  and  Pladbtm  Core  Energy  Versus  111116 


Hgure  (34)  Total  Energy  and  Flatform  Core  Energy  -  First  100  Seconds 
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Case  5:  bAt<l,  Stable,  Verfaul^  Model 
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gufe(57)  Modded  and  Actual  Flatfonn  Core  Energy  and  Nutation  Angle  Versus  Tu 


Figure  (58)  Modded,  Actual  Platform  Core  Eneigy  and  Nutatkm  Angle-  Rrst  100  Seconds 
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Figure  (S9)  Core  Enogy  Versus  Nutation  Angle 
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CASE  6:  ^  <  1  ,  Unstable,  Verhulst  Model 


Systeoi  PsTuncters 


Platfionn 

Rotor 

V-  1600  kgm2 

/i- 

1000  kgm2 

1500  kgm2 

/3* 

1200  kgm2 

hr»  1000.0  kg 

700.0  kg 

1.0  kg 

m  = 

1.0  kg 

cf  »  1.0  m 

as 

1.0  m 

A'-  1.0  S 

m 

As 

m 

1.0 

sec 

Cs 

1.0  JS- 

sec 

L»  1.0  m 

htoui 

*  0.897 

Initial  Conditioas 
Platfbnn  Rcxcr 


/  s 

0.0  m 

Z  s 

0.0  m 

i' s 

0-0? 

i  ■ 

0.0  m 

Oil  • 

7.27  X  10^  nd 
sec 

<»»» 

1.5 

sec 

a>i  s 

0.100*1^ 

sec 

0)2s 

0.0 

sec 

Gne  Energy  Panuneters 
Platform  Rotor 


Ec  initial  ~ 

15.1  J 

Ecinitial’^ 

3061.6  J 

Ecfmai  = 

.1161.0  J 

Ecfmai  * 

1489.9  J 

A' 

positve 

&  = 

A 

p(»itive 

r'  s 

-1x10-®  s-i 

r  = 
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Figure  (60)  Total  Enogy  and  Peroentl^erenoe  Angular  Momentum  Versus  Time 
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Figure  (61)  Total  Energy  and  Fcroent  Difference  Ai^ular  Momentum -Hrst  100  Seconds 
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Figure  (62)  Total  Energy  and  Rotcr  Core  Eneigy  Venus  Time 
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IQS 


concocfiy  (J)  I  coweoeriy  (J) 


enerfy(J)  I  eiMr|y(J) 


coreeoergy  (J)  I  coieeiwify  (J) 


cofeeiieify(J) 


Cue  6:  b^Kl,  UnsUble.  Verhuist  Model 


0  1000  2000  3000  4000  5000  6000 

tiine(Mcaods) 


Figure  (66)  Modeledand  Actual  Rotor  Core  Energy  and  Nutation  Ai^e  Versus  Time 
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Figure  (67)  Modded,  Actual  Rotor  Core  Energy  and  Nutation  Angle  -  First  100  Seconds 
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ngure(73)  Total  Enogy  and  Rolor  Cere  Enogy  Versus  Tune 


Rgiire(74)  Total  Energy  and  Rotor  Core  Eneigy- First  100  Seconds 
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Figure  (76)  Total  Energy  and  Platform  Core  Enagy- First  100  Seconds 
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Figure  (80)  Modded,  Actual  Flatform  Gore  Enogy  and  Nutation  Angle-  first  100  Seconds 
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APPENDIX  B  COMPUTER  PROGRAM  CODE 


/ 


declarations 


/ 


•ineluds  <Mth.h> 
fineluda  <nMlloc.h> 
lineloda  <stdie.h> 

•daflna  MMSTP  ISOCSO 
fdaUna  TINY  1.0c><C 


/•  odcin'  '/ 
/•  o<lpiii<  • 


tdatina  PCROH  -0.20  /*  •/ 

•dd£ina  fSIWNK  -0.25  { 

«da<ina  rCX>R  0.0066«666<6<CC06666C666<666(66667  /*  l.O/lS.O  keqe*/ 
«da£ina  SAFETY  0.9  /*  •''V 

•dafina  EMtCOK  6.0a-<  /*  '/ 


fXdat  **y*0,  *xx«0:  /*  dafinlng  declaration  rkduaih  */ 
int  kmax«0,  kount-0;  /*  daCining  declaration  odeint  ♦/ 
£loat  •xp»0,  ••yp-0,  dxsavO;  /*  defining  declaration  odeini  •/ 


/••••••••••*•••••*•  error  function  •••*••••••••»*•••••••*/ 

void  nceccoc(ercor_textl 
char  error  text  I ] ; 

1 

void  exit ( ) ; 

fpsintf (stderr, "Numerical  Recipes  run-time  error ... Xn”) ; 
fprintf  (stderr.  "%s\r.‘*, error  text)  ; 
fprintf  (stderr.  **. .  .row  exiting  to  system. .  .Xn")  : 
exit (1) ; 

) 

/••••••••»**********c«*  vector  function  »•»*••*••»*»••••*•******/ 

float  *vectorr (nl,n)t) 
int  nl,  nh; 

{ 

float  *v; 

v»  (float  *)maiioc I  (unsigned)  (nh-nl-^l)  *siseof  (float) ) ; 
if  (iv)  nxerzox (“allccation  failure  in  vector ()”); 
return  v-nl; 

) 

/••••••••••••*••#•*«•  matrix  function  . . . 

float  **matsix (r.zl. r.r'r.,  ncl, nch) 
int  nxl,  nx)i,  nci,  .*<rh; 

I 

int  i; 
float  ••.ti; 

m  ■  (float  ••)  malic: ( (unsigned)  (nr)>-nrl+l) *sixeof (float *)) ; 
if  (!m)  nxerrox ("allocation  failure  1  in  matrix ()"); 
m  —  nrl; 

for  (i«nxl;i<-nxh;i'»-) 

( 

m(i]  •>  (float  *)  ma Hoc  ( (unsigned)  (nch-ncl'*-!)  *sixeof  (rio.-ii ) )  .- 
if  (imti'j)  nrcxroz ("alloc  tion  failure  2  in  matrixO”); 
m(i)  —  ncl; 

) 

return  m; 

i 

/•••••••••••••*••.•••  «jee  vector  function 

void  frec_vector (V. nl.nh) 
float  "v;"" 
int  nl,  nh; 
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» 


/•*••••******••••••••*•••  CxM  Mtrlii  function  . . . 

void  fc««j«otclx(m,ncl,nch,nel,nch) 
float  **«; 

lot  nxX>  nzh,  nel,  nch: 
t 

int  1: 

£or(l-nrh;i>»ncl;i— )  froaCCchac*)  (mUHncD); 
fxo«(  (ehas*)  (in^nrl) )  ; 

\ 

/••*••••••*••••*••••*••  xk4  function  ***•*•*•**•*•*«**«*«♦•/ 

void  ck4ty>dydx,n.x.h,yout,darivs) 

float  ytl,  dydxd,  x,  h.  youtd: 
void  (*dorivs) () ; 
int  n; 

int  i: 

float  xh,  hh,  h6,  ‘dyin,  *dyn,  *dyt,  *yt,  *vectorr(); 
void  f zo*_voctor ( ) ; 

dym  «  voctozr (l,n) ; 
dyn  •  v«etocc(l,n) ; 
dyt  ••  voctocr(l,n); 
yt  -  voctotr (1, n) ; 
hh  -  h*0.9; 
h<  -  h/6.C; 
xh  ■  x-»hh; 

for  {i-l;l<-n;i+4)  yt(i)  -  y (i]4hh*dydx|i] :  /•  first  stsp  •/ 

(*doriva) (xh,yt,dyt,dydx) ; 

for  (i-l;i<-n;i44)  yttil  -  y{il4hh«dytlil; 

(*doriv«)  (xh,yt,dyir.,dyt) ; 
for  (i*l;i<»n;i44) 

( 

yt(ij  -  y(i)*h»dym(ij; 
dyn(i]  ••  dyt(i]4dym(i]; 

(•dtrivi) <x4h,yt,dyt.dym) ; 

for  (i-J;i<-n;i44)  youtli]  -  y (i]4h6« (dydx(il4dyt (ij 42 .0 ‘dyn l i i )  . 

froo_voctor(yt,l,n); 

ftoo^waetor (dyt . 1 , n)  ; 

f  *oo_v*ctor  (dym,  1 ,  n)  ; 

frao^vtetor (dyn, X, n) ; 

) 

f******************’***  rkdumb  function 

void  rkdumb (vstart, nvar, xl, x2, nstop.derivs) 
int  nvar,nstop; 
float  vttartIJ,xl,x2; 
void  (*d«rivs)  0 ; 

( 

int  i,  k; 
float  X,  h; 

float  *v,  ‘vout,  ‘dv,  *voctorr(); 
void  rk4 ( ) ,  nr error ( ) ,  f r o«_voctor  < ) ; 

V  -  voctorr (l,nvar) ; 
vout  »  voctorr (l,nvar) ; 
dv  •  voctorr (l,nvar); 
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Cox  (l"l;i<»nv«x,*i^*) 

( 

vCU  -  vstxrtlH; 
yliJCD  -  vlil; 

I 

xxlll  -  xl; 

X  -  xl; 

h  *  (x2-xl)/nsCxp: 

Cox  (k«l:k<«nst«p;k'»'») 

( 

(*d«xivs) (x,v,dv) ; 

xk4 (v,dv,nv«c,x,h,vottt.d«xlvs) ; 

If  (x^h  “  X)  nxcxxox ("Step  site  too  axMll  in  toutin'*  nKi'i’nn") 

X  h; 

xxlk'*'!)  -  x; 

for  (l»l;i<-nw«r;l+*) 

( 

v(i)  •  vout (ij; 
ylijlk*ll  -  v(ll; 

) 

» 

fxoo_v«ctor (dv, 1, nv«r) ; 
f xoo~voctox (vout, 1, nvxr) ; 
f roo~voctor (v, 1, nvar) ; 

) 

/•••••••••••••••••**••  rkqe  function  **••••*•*•**••***•••*•«••*••/ 

void  xkqe (y, dydx, n, x, htxy, opx, yscal, hdld, hnoxt, dot! vs) 

float  yd,  dydxd,  *x,  htry,  aps,  yscald,  "hdid,  *hnext; 
vold  ( "darlvs ) ( ) ; 

Int  n; 

( 

Int  1; 

float  xaav,  hh,  h,  tamp,  aremax; 
float  *dyaav,  *ysav,  *ytamp,  "vaetosxO: 
void  ck4  ( ) ,  nraxxot ( ) ,  f xaa_vaetoc ( ) ; 

dyaav  •  vactorr(l,n); 
ysav  •  vactorr (l,n) ; 
ytanp  •  vactorr (l,n) ; 
xsav  -  {‘x) ; 

for  (i-l;l<-n;i**) 

( 

yaavlij  -  y(ij; 
dysav(i)  -  dydx(l]; 

) 

h  <■  htry; 
for  (  ;  ;  ) 

{ 

hh  -  0.5'h; 

rk4 (ysav, dysav,  n,  xsav, hh, ytan^, darivs) ; 

•x  -  xssv'fhh; 

(*darivs)  (*x,ytainp, dydx); 

rk4  (ytaiiq),dydx, n,  *x,hh,y, darivs); 

•x  •  xsav4h; 

if  (*x  ••  xsav)  nraxxor  ("Stap  sisa  too  saiall  in  routine  RKOC:'') 

xk4  (ysav,  dysav,  n,  xsav,  h,  ytaxf),  darivs) ; 

arxmax  •  0.0; 

lot  (l-l.-K-n;!*^) 

( 

ytampli]  -  y(l]-yteiiv(i); 

tamp  >  fobs (ytampli] /yscal(i)): 

if  (arrmax  <  tamp)  arrmax  •  tamp; 
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•cnwM  /«  ffps; 

1C  (•rzMx  <*  1.0) 

( 

•hdld  •  h; 

*hn*xt  -  (•rXMM  >  CRRCON  7  SKFeTY‘h‘ftxp(PGnow*  (•mj  (oi  im  o  M 

1  .  (1  • 

break; 

) 

h  -  SM‘ETy*h*«Mp(PSIUWK*l09(«rtMx)); 

Cor  ;'!♦♦)  yfi)  ♦-  yt««pCi)*FCOR: 

Croo_voetox  (ytoiap,  1.  n) ; 

Cro«_v«ctor (dysov.  1 ,  n) ; 

Cxoo_v«ctor (yoav.  1.  n)  ; 

»_ 

/•••••••••••*•••••.•«  function  •««••••*••••»»•«•••••/ 

void  odoint  (y start .  nvar , xl,  x2 ,  ops ,  hi ,  hmln,  nok,  nbad,  d« c i us ,  t kci<- ) 

Cloat  ystartf],  xl,  x2,  ops,  hi,  hmin; 
int  nvar,  *nok,  •nbad; 
void<*dorivs) () ; 
void(*r)(qe)  () ; 


( 

int  nstp,  i; 

float  xsav,  X,  hnaxt,  hdid,  h; 

Cloat  •yscal,  'y,  *dydx,  •voetozrO; 
void  nrarrorO,  f raa_voetor  ()  ; 


yscal  -  voctorrd.nvar); 
y  •  voctorr (l,nvar) ; 
dydx  -  voctorr (1, nvar); 

X  >  xl; 

h  -  (x2  >  xl)  7  Cabs (hi)  :  -Cabs (hi); 

•nok  -  (*nbad)  »  kount  -  0; 

Cor  (i-l;i<-nvar;l4*)  y(ij  -  ystartti); 
iC  (knax  >  0)  xsav  •  x>dxsav«2.0; 

Cor  (nstp-1 ; nstpoKAXSTP;  nstp^-k) 

( 

(•dorivs) (x.y,dydx); 

for  (i«l;i<-nvar;i-»-»)  yacal(i)  -  tabs (y (i) ) +fabs  (dydx I  i )  •!.)  •  TINY, 
if  (kmax  >  Cl 
( 

if  (fabs(x-xsav)  >  fabs(dxsav)) 
i 

if  (kount  <  kmax-l) 

( 

xp(4"«’kount)  "  x; 

Cor  (i«l;i<-nvar ;!•'*')  yp(i)  (kount)—/ U)  : 

xsav  •  x; 

) 

) 

) 

if  ((x-»h-x2lMxfh-xl)  >0.0)  h  -  x2-x; 

nvar ,  tx,  h,  ops,  yscal,  Chdid,  Chnoxt ,  dcr ivs)  ; 
if  (hdid  —  .*:)  -M-(*nok);  also  •’•(•nbad); 
if  ( (x-x2)^ (x2-xl)  >>0.0) 

( 


for 

if 


(i>l;i<>nvar;i4>) 

ystart(i)  -  y{ij; 
(kaiax) 

( 


xp(*’»kountJ-x; 

for  (i>l;i<>nvar;i^^) 


yp(i] (kount) 


y  t  i  I  ; 
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I 

(r«t_v«etor (dydx, 1, nwar ) ; 
fr«t~v«etor (y, 1, nvar) ; 
f x»t2««etor (yseal . 1 , nvar ) ; 
taturn; 

) 

if  (faba(hnaxt)  <*  loalii)  ntazror ("Stap  alsa  too  siull  in  oukiiit’*). 
h  •  hnaxt; 

) 

nzaczor ("Too  aiany  stapa  in  zoutina  OOCINT") ; 

) 
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. . 

•iiMltid*  ”Klik.li* 


i 


•daXiM  M  7 

«4aflM  NMCMWAY  100000 
fdaXlM  SOU)  (U)*U)) 

lat  l(0iint,  kmm*t 
float  **yp«  *«p> 
float  dxaav.  daain,  dxlliait; 

float  Zl-0.0,  Z2-0.0.  Z9-0.0.  M-0.0,  at-O.O,  nJo^O.O,  a>0.0,  k-O.R,  r-n.;); 
*l®*t  ZlpH).0, Z2p*0,0, Z3p*0.0«Hp«0.0,ap«0.0,aibp«0.0«ap-0. 0, kp-O . o . •»; 
float  1«»0.0,  slana  0.0«  MT«0.0,  v»0.0,  Kho*0.0,  rltop*0.0,  Ll-o.o, 
float  A-0.0,  C-0.0,  JJp-O.O,  w3p-0.0; 

. . •••*•*.•  CQOATZONS  or  HOTION  . . . . . 

void  dociva(K,y,dydx) 
float  *.  yll,  dydxt); 


float  sata»0.0,  dydxxata*0.0; 


float 

float 

float 

float 

Al-0.0, 

A9-0.0. 

A17-0.0, 

A2S-0.0. 

A2-0.0. 
A1 0-0.0, 
AlS-0.0, 
A2  C-0.0. 

A3-0.0, 

All-O.O, 

A19-0.0, 

A27-0.0, 

A4-0.0, 

A12-0.0, 

A20-0.0. 

A2S-0.0: 

A5-0.0, 

A13-0.0, 

A21-0.0, 

A6-0.0, 

A14-0.0. 

A22-0.0, 

A7-P.0, 

A15-0.0. 

A23-'0.0, 

Aft-n.n 
AJfi-O.O 
Ari-n , 0 

float 

float 

float 

float 

Bl-0.0, 

B9-0.0, 

B17<-0.0. 

B2S-Q.0. 

B2-0.0, 

BlO-0.0, 

BlS-0.0, 

B2C-0.0: 

B3-0.0, 

Bll-0.0, 

•19-0.0, 

B4-0.0, 

B12-0.0, 

B20-0.0, 

BS-0.0. 

B13-0.0, 

B21-0.0, 

B6-0.0, 

B14-0.0, 

B22-0.0, 

B7-0.0. 

B15-0.0, 

P23-0.0, 

PB-O.O. 
ni6-o.n.- 
K," ;  n.O; 

float 

float 

Cl-0.0, 

C9-0.0, 

C2-C.0, 

ClO-0.0, 

C3-0.0. 

Cll-0.0; 

C4-0.0. 

C5-0.0, 

CC-0.0, 

C7-0.0, 

c^n.it; 

float 

01-0.0. 

02-0.0, 

03-0.0, 

04-0.0, 

OS-0 .0, 

DC-O.O, 

07-0.0, 

PR -0.0; 

float  Cl-O.O, 
float  C9-0.0, 

E2-C.0. 

SlO-0.0; 

E3-0.0, 

E4-0.0, 

E5-0.0, 

E6-0.0, 

E7-0.0. 

KR-n.O; 

float 

tl-0.0. 

t2-0.0. 

13-0.0, 

E4-0.0. 

ri-0.0. 

r2-0.0; 

dydxUJ  -  yC5J; 


dydxIO]  -  y|7]; 


AX  - 
A3  - 
AS  - 
A7  • 
A9  - 
All  - 
A13  - 
AIS  - 
A17  - 
A19  • 
A21  • 
A23  - 
A25  - 
A2C  - 
A2t  - 


-  sho*y(4]«zhop*y(6]; 

-  -(A-C)*y(2J*y(3J; 

■  2*nT*sata*dydxzata*y(l] ; 

•  -HT*SQ(sata)*y(2)*y(3]; 

■  2*m*(«ata<»1.2)*yM)*y(2]*y(3]; 

«  -2*m*y(4)*dydxt«ta*y(l]; 

•  ■•SQ(yMJ); 

■  >2*ap*(sata-Ll)«y(€]; 

•  •2*ap*(zata-U)*y(7]*y(l); 

■  2*»P*yI«l*yl7J*y|ij j 

•  •»p*S0(y(C))*yi2]*y(3]; 


dydxsata  -  rho'ylSJ-iihpp*-,;’); 


A2  - 
A4  - 
A(  - 
AS  • 
AlO  • 
A12  • 
A14  • 
A16  • 
AIS  • 
A20  - 
A22  - 


J3p*si9ai«*y(2}; 

MT*SQ((ata) ; 

-2»m»  (saca«L2)  ‘yMJ  ; 
-2*iD*U«ta4^L2)  ‘ylSJ  ‘yll )  : 
2*m*y(41*y(51*y(I]; 
-«*SO(y|4))*y|2J‘y|3J. 
-«*yt41*a*ylll»yl2l; 
2*np*(s«ta-Ui  *y|6]  •y|2J‘i  (3I.- 
~2*ap*y|6)*dydxiotn*y 1 1 j ; 
i<V*SQ<yl6)); 

“•p*y|6J*ap»co8(siqTi*K)  ; 

A24  “  xip'ap*!; > !i (s inn'-')*x) ; 


IT.  -  -xa»-y ivj -ap-cos (SIM 

-^•y]SJ*ap*eoa<tl9aw*x)*yll)*yl2);  A24  -  xVapT.i 

^*ap«ain  (aigoa'x)  -y I S)  •  (SQ «y  1 3)4si9iaa)  -S0(y  12} ) )  ; 

?T^«**t*'i*^^*‘^**®*  *  A13+A22; 

A14A24A34A54A74AS4A94A104A124A144A1€4A174A184A194A214A23^A7'; 


SI  -  -(C-A>*yfiJ*yl31; 

B3  •  2*KT*aata*dydx<tta*y(2] ; 


B2  -  -J3p*al9fna*yl  1]  ; 
B4  -  MT*SQ(xela); 
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B5 

>7 

■9 

Sll 

B13 

•IS 

•17 

•19 

•21 

•22 

•24 

•2C 


MT*SQ(t«t«)*y(ll*yC3]:  •• 

-2*»i*  (««t«4L2 J*yl41»ym*y(31;  •• 

-3*»*y (4) *dydxi«ca*y (2] ;  BIO 

«*SOIyl4});  Bll 

-■•a;  ^14 

-2*ap*(tata-Ll)*y(*J;  •!« 

-2*iiip*(»ata-l.H*yni*yl2);  •!• 

av*SQ(y(C]);  B20 

-iiip*ap*cot  (sl9iu*x) ;  B23 


-2*«i*  (*ata«1.2)  *y  M  !  ■' 

-2»i«*  data*  1.2) 
2*ai*y|4)*y|51*y|7I; 
Bi*SO(y|4))»yll)‘yl'»,- 
-ai*a*y|41*(SO(yl31 1  -rO(vM  I ) )  • 
>2*iip*  (tata-L) )  ‘y  1 ‘ y  m  •  V I  ■<  I . 

2*aip*yl«|*y)21*  )yn  -.iy<).;.  >M  .,) : 

BV*SQ(yl61) ‘yUPy  I 
-ap*ap*ain(si9a<'** X)  *)  ifil; 


••p*ap*cea(si9BM*x)  •(SQ(y(3)‘*'Sl9M)-SQ(y(l)) )  *y  (6|.- 

np*ap*alh(si9iaa*x)*y|C)*yilj*y|2);  B25  -  A*-n4*n0‘ =  1 )  mm  v-m 

•l-*B2^B3-*BS-*B7«B84B9^B104B12«B14*-Bltm7^B19«B20«B22*D24; 


Cl 

C3 

cy 

C7 

C9 

CIO 


-2*m*a«yI51*ym;  C2  -  -■••a^yMJ; 

«*a*yl41*y(21*y(31;  C4  •  -2»ap»ap»aiB(tl9m8‘:-.)  *y  I  M ‘yl’l ; 

-aip*ap*aiii(ai9Ma*x)*y(C];  C«  -  -aip*ap*aln(aigiaa*x)  •>- 1 r) ‘y  1 1 1 ‘y  I  '1; 

-2*Bip*ap*eos(si9BW*x) *y(7) *yll);  Ct  -  -■p*ap*coaC8iq»a*x) 'y K>l  ; 

ii^*ap«eoa  (sigaM^x)  ‘y  ( 6]  *y  121  ‘y  (3|  ; 

C2>C8;  Cll  -  C1+C3-*C4*C6*C7*C9; 


D1 

D3 

05 

07 


m«(l-tho);  02  -  -i»p*t)»o; 

-m*a;  04  »  ai*a*y(l)*yI3] ; 

-m*  (SO  <y  ( 1 1 )  *80  (y  12 1 ) )  *  (y  1 4 1  *  (1-the)  -t2-rl»op*y  ( 6 ))  ; 
)c*y(4];  08  -  D4*^05*0C-*D7; 


O'-  -  r-ylM: 


£1 

£3 

£S 

£7 

£9 


-m«rhop;  £2  •  aip*  (l-r)>op)  ; 

iap*ap*ain(alg>na*x);  £4  •  aip*ap*ainC8igma*x)  *y  12)  *  (y  3)  •  ?' " i'.n"  M 

-inp*ap*eoa  (sigma'x) ;  £8  -  ap*ap*coa(sigma*x)  *y  ID*  ty !  Sl+J'r.ioir  m 

-nip*(SQ(y(l))-*SQ(yI2)))*(y(6)*(l-rhop)-*Ll-rl»o*y|4));  TP  -  rr'ylM 

kP'ylS);  £10  -  E4+£6*-E7+E8+E9; 


21 

22 

ri 

23 

24 
£2 


(- (A27*B13) / (A28*B23) >£l/£3) / ( (A27*B2S) / (A2S*B23) ♦E5/E3) ; 

(•  (A27*B21)  /  (A2C*B23)-*A24/A2«-C2/£3)  /  ( (A27*B2S)  /  (A26*B23)  *F.S/T.3) ; 

(- (A27*B26) / (A28*B23) ♦A20/A26-£10/£3) / ( {A27*B25) / (A26»B23) *  E3/E3) . 

( (C*B13) / (C10*B23) • (A27*B13) / (A2C*B23) ) / 

(  (A27«B2S) / (A28*B23) - (C*B25) / (CIO  *  B2  j ) *CS/C 10 ) : 
( (C*B21) / (C10*B23) - (A27*B21) / (A2C*B23)^A24/A26) / 

((A27*B25)/|A26*B23)-(C*825)/(C10*n23)*Ci/C]C) . 
( (C*B26) / (C10*B23)-Cll/C10- (A27*B26) / (A26*B23) ♦A2e/A26) / 

( (A27*B25)  /  (A28*B23)  -  (C*B25)  /  (CIO *1’/' 3 )  '  Cf'/C  10). 


dydxlS]  -  ( (-F2-D8/D3)/(24*-02/D3)4(Fla08/03)/(Z2-*D2/D3) )/ 

( (-21-01/03)  /  (2202/03)  ♦  (23+01  /P3 )  /  I  r  1 1 1>?  /»))); 

dydx(7J  -  ((-r2-D8/03)/(23+-01/D3)+(Fl«08/03)/(21+01/03))/ 

( (-22-02/03) / (21+01/03) +  (21*02/03)  /  C  SMII  IM)  ); 

dydx(31  -  (dydx(S]*(CS«B13)/(C*B2S)+ 

dydx 1 7 J • ( (C10*A24 ) / (C+A26)  +  (C5+B21 ) / (C+B25) )  + 

(C10+A28)  /  i;C*A26)  ♦  (CS*B2«)  /  (C+B25)  -Cll/C)  / 

(1- (C10*A27) / (C+A28)- (C5*B23) / (C*B25) ) ; 

dydxlD  -  dydxt3]*(-A27/A26)+dydx|7]*(-A24/A26)+(-A28/A26); 

dydxI2]  -  dydx(3)*(-B23/B25)+dydx(S]*(-813/B25)+dydx(7)  *(-B21/D?:<)-:-«2ri/l«;  r>; 


ppoGrah  . . . 

nainO 

I 

Int  i.  ),  nbad,  no)i; 

int  11-1000,  12-2000,  13-3000,  14-4000,  lS-5000,  16-6000.  17-7000. 

Int  18-8000,  19-9000,  110-10000; 

float  aps,  M,  (min,  dxsavd,  xl,  x2,  xacala,  *ystart; 

float  Hl-0.0,  H2-0.0,  H3-0.0,  fllp-0.0,  N2p-0.0,  H3p-0.0.  tcm-O.O.  gum'  v.- . 
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eiMt  TlHI.e.  T2-0.0.  T9-0.0.  T5-0.0.  T«-0.0,  Tf~0.O, 

noat  Tt«4.0.  TXO-0.0,  Tll-0.0.  T12-0.0.  T13-0.0; 

CXoat  sflKifld.  Xt,  Xtp,  Xttatal.  la,  lap,  latotal,  wn,  lambd-t, 

XiMt  *ka,  *kadal,  •atetal,  •atotaldal,  *h,  •hdal,  ‘titata; 

float  •Kc,  •tedal,  •Cep,  *£^681,  •Cehypa,  *Cchypadel.  *Ccphyp<',  'r<piiv|v.|. 

float  0,  Op,  *eta,  •atap,  Ceadn,  eefsctor,  acpfactor,  haq,  Ecfinii,  c.-pi in 

float  aipnoe.  ai«aoep,  ai«nq,  slqnqp,  at«,  arqp,  Eef,  Ecfp; 

float  •Eeoap,  •Cepaxp,  •Eeaiipdol,  *EepoxpOal,  •atah,  *fftahp; 


ystact  -  vaeterr (1,M); 
ka  -  vaote»(l,NM(M(MY).< 

kadal  -  vaototz (l.MAXACCAy) ; 
atotal  -  aaeterzd.HAXMUMy); 
atotaldal  -  aaetotrd.MkEAMlAY) ; 
k  -  «aeterzd,MM(MUtky).- 

thata  -  aaecotcd,HM(MUMY); 
hdal  '  -  aaetotrd.MMCMAAY); 
Ceaxp  *  vaetetzd,HaXkmAY): 

Cepaxp  -  vaeterrd,MO(AmAY); 


Ec  “  vector r ( 1 , "AXAruAV ) 

Ecpdal  >  vector  r (1 . maxarhay ) 
Echypa  •  voctorr  (] ,  tiA.XARKAY) 
Echypedal  •  vectorr ( ] .  ttAXAPRAVi 
Ecphypa  •  vectorr  (l,:!''XAUi»AV) 
EcphypadaX  >  vectorr  d.nr.X.VRRAY) 
ata  -  vectorr (l.rjAXARRAY) 

atap  •  vectorr  (l.r:;.XAP.RAY) 

Ecaxpdal  -  vector r  ( 1 ,  tt;- xabray  ) 
Ecpaxpdel  -  vectorr (1,;;;.xakpay) 
etah  -  vectorr (l,;^^:■:AKRAV) 

etahp  •  vector  r  ( 1 ,  rv.XARRA y ) 


acanf <"%f%f%f%f%f%f%f%f%f%f%f%fif%{%{%{%f«(tf%ftctftftf%fiftfT'* 

ax2, (XI, (13, (M, (o, (a, ck,(c, (L,(Ilp, C13p, tMp, top, cap, ckp, ^ cp. 
4**3p,  (ystazt  (1 J ,  (ystact  C2 1 ,  (ystart  J 3 J .  Cystart  H ) .  ty st a r •.  i  5 ) , 
•yakact C<J < (ystact (7] , Saps, (acfactoc, (acpfactor.  tEcf ,  (E-rp) 

xl-0.0001;  hl«0. 00001;  haiin-l.Oe-11;  kMax>100000: 

*b  -  m;  obp  .  mp:  12  -  II;  I2p  -  ixp; 

•i9M  •  a3p“ystart{31; 


/•**•  critical  daapinp  -  input  danpinq  -  XOOO.O,  makes  it  crit  -M  •*••/ 
If  (e  «-  1000.0)  e  -  X.3*S4ct((*ia*k); 
if  (ep  ••  XOOO.O)  ep  -  l.S*sqct(4*mp*kp); 

/•••  establish  tims  interval  for  savinq  data  *••/ 
dxmin  ••  1.0a>(; 
dxlimit  -  (x2>xl)/17S0.0; 
dxsav  -  (x2<20.0)  7  dxmin  :  dxlimit; 


HT 

V 

A 

C 

J3p 

z))o 

chop 

LX 

L2 


M^Mp*4 .0*ml>«4 .0*Nbp; 

(Hp«4.0*i*p)  /MT; 

Xl-kXXp^X . 0*Rb*SQ (a)  +2 . 0*mbp*SQ (sp)  *  iHp*4 . 0 ‘mbp)  *  U  •  0-v ) 
I3>I3p-(4 . 0*fflb*SQ  (a )  *4 . 0*iiibp*SQ  (ap)  ; 

I3p+4 . 0*aibp*SQ  (ap)  ; 
m/MT; 


mp/MT; 

L«  (H-(4.0*mb)/MT; 

L«  {Mp-»4 . 0*mbp)  /MT; 


i:.)  ; 


ccmriqid  -  ((Hp'f4.0*mbp)*L)/MT; 

It  -  XX'fH*SQ(<emciqid)'»iiib*(4.0*SQ(ccfflciqid)4'2.0*SO(a) )  ; 

Xtp  -  XXpaiq>*SQ(L-scmriqid)-»inbp*(4.0*SQ(L-scmciqid)+2.0*SC'-  ) ) : 

IttotaX  - 

Xs  -  X344.0*R<b*SQ(a); 

Xsp  -  X3p>4.0*mbp*SQ(ap) ; 

Xstotal  -  Xsalsp; 


®4^^at (ystart, N, xl, x2,aps, hi, hmin, (nok, (nbad,  dec i vs,  rkqc) ; 
for  (j-X;  j<>kount;  j4"f) 

tern  -  ((Hp^4.0*mbp)*L«m*yp|41 {j]*np*yp(S) ( j))/MT; 

HX  -  (21«H*SQ(scm)-Ha*(2.0*SQ(a)«4.0*SO(scm)4SQ(yp|4) I > 
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■2 


•2.o*tai*ypi4]  I  ji))*ypUl  I j)-»*«*ypi<)  I H 'ypi^i  I  Jl* 

-  (X2-«l»SQ(iai)'Hi*(2.0*80(«J*«.0*SQ(«e*>*SO{ypMn jn 

-2.0*tai*ypl«JIH>)*ypl2H  j)-»*a*ypl51 1  Jl  ■ 

■3  -  -••••ypl41 131*yplll  ljJ  +  (I3*«.0*»*SO(«))*ypt3J  IJl-- 

■Ip  -  (np*HP*»0(*>-«oa)“Hi**(2.0*SQI»p)M.0*80(t-«c»MSQ(ypj6J  I  jj)< 

2 .0* (t-soi) ‘ypl*!  C31 » > •yptlHll-^*»P*coa jn  • 

»pj«j  I  ji-ypcij  I 

'  yp(C]  (j]<H9*ap*«lii(si9M*vpi  j))*ypn)  I  j).* 

I2p  -  (X2p*np*SO(L-i«i)-Hip*(2.0*SOI«P)^<*0*SQ(t-te*)*SO(ypl6)  I  jn  ' 
2.0*(X*-«ea)*ypl<]  I  jl))*ypl2J  I  jl-«»*»P*«in(»i?«*>«Pl  Jl » ‘ 
ypl«]  1 1) •yp(3]  1 31  -in»*«p*«igii«*»in  U » * 

ypIC]  I  ji-Hp*ap*eea(ai9M*xp|  j])*ypp]  |j]; 

■3p  - -■p**p*eo*<*i»“**pi3n*ypi*HJJ*ypi^i •  ... 

•lB(alfl**iep|jl)*ypl«](3J*ypC2JI31*II3p*«.0*i«^>‘SO<*p))* 

(yp(31(j]'»«i«M)-‘ 

htjl  -  •qrt(S0(Hl4Hlp)>S0(B2-HI2pl4S0(a3^«3p)>; 


•rguMiit  •  (H3-»M3p)/h|  j]: 
if  (ar9unant>0 .99999999)  arguMiit  •  1>0; 
th«ta(j]  •  S7.2957795131*acos (argufliant) ; 


hdaKj]  -  h|3J-h(ll; 

n  -  0.5*(Xl«SQ(yptinjl)4X2*SQ(ypl2)l31)+I3»SQ<yp|31tjJ)); 

T2  -  0.5*(Xlp*SQ(yp|injn*X2p*SQ(yp|21ljJ)*I3p* 

S0(yp(3)  I  J)<r.>'jmK  )  .■ 


t3  -  0.5*ia*lSQ(ypl5H31)^(S0<ypll)lj))*S0(ypl2njJ))‘S0(ypl'lll.il> 

•2.0*a*yp(l) t jJ*ypI3) I jl ‘ypli ) I j) ) ; 

T4  -  0.5*BV*(SQ{ypI7](  j])*SQ(ypt61ljnMSQ(ypll)  I  jl)  ♦S0(ypl21 1  j) ) ) 
♦2.0*ap*ain(aigma*xpl j] ) *ypl2) I jl •yp|3) I i) *ypl6 j | j) 
-2.0*ap*eoa  (aigma*xpi  j]i*ypil  j  I  jl  •ypl3)  I  j)  *yp|5)  I  jl ) 

T5  -  -ia*a*yp[2Hjl*yp(5njl; 

T6  -  -a4>«{-8p*ypI7n j)«yp(21ljl*coa<aigma*xp|j)) 

•♦ap*ypi7) ( j]*yp|li j ji*ain(aigma*xpl jl) 

•ap*algBia*yp(ll ( jl*ypl61 I jl*coa(aigma*xpl j)) 
>ap*aigma*yp(2i j jj •ypi61 j jl*ain(aigma*xp| JJ ) ) ; 

T7  -  (0.5*SQ(m«yp(51Ijl*i«p*ypl71(jl))/MT; 

T$  -  (0.5*(SO(ypllHjl)*SQ(ypl21Ijl)) 

•SQ  Cia*yp  1 4  H  j  1 -•  mp*  yp  1 6 1 1  j ) ) ) /in  ; 

19  -  0.5*(SQ(ypIll  (jl)+SQ(ypl21l31))*l»«>4.0*a»)»SO(L)« 

SQ(  (Mp4'1.0*RT))/MT)  ; 

TIO  -  0.5«(SQ(ypIll  [jl)+80(yp(21C31))*<»<P*<-®*»P)*SQ(X.)* 

SQ((M+4.0*in)/MT); 


Til  -  <iii4ag»*(-ypI511jl)Mai*ypl51tjl+ag>*ypI71|j))/MT; 

T12  -  -m*yp(41Ijl*(SQ(yp|ll[jl)*SO(ypl21Cjl))* 

(ia*yp(4)  I  jl*aip*ypl61 1  jl-»  (Hp+4 .0*nip)  M.)  /HT; 
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m  • 


-*i»*ypl«)lJI*<»0Cy»UlIlll*SQIyp|2HJin* 

(■•ypMI  ()H«p*yp|fi|  I  •i.)  /ni 

k«IJI  •  T1*T2*T34T<*T5*T«*T7»T8<T»*T10+T114T12'T1  1; 

•tetaKlI  -  k«tJl*0.5*k«SQ(ypl<Mjn*0.5*kp*SQ(ypK.|  I  ill ; 

•totaldvll)]  «  atotall j] '•total ID : 

Eelj]  -  O.S*(Ittotal*(SO<yplDlllDSO(ypl2Jljl)) 

♦  Iatotal*SO(ipni  I  il 1 1 

Cedolljl  •  CcIji'CelD; 

Ccpij)  -  0.5*(Xttotal*(SQ(yplDljlDS0(yp|2]ljl») 

'»lstotal*SO(ypl 7)  I  il  iriamii  i 

Xcpdallj]  -  Ecpljl-Ccpll].* 

» 

wn  ■■  (Xa*yp{3)  {l]-»Itp*<ypt3)  (D-^siqiwD/Xttotal; 

laadada  -  «m-ypl3]|D; 

laaMap  -  «m-(yp(31  lU^al^M)  ; 

baq  -  SQ(Xttetal)*(SO(yp(l]  (in-»SQ(yp(2nDD 

♦SQIls^ypO)  UDIap*  (ypl3)  I '  I  'r.iom-*) ) 


iC  <th«ta(keuntl-th«ta(l]  <  0.0) 
I 

if  (lambda  <  0.0) 

•la« 


•laa 


if  (lambdap  <  0.0) 

•laa 


if  (lambda  <0.01 

•!•• 

if  (lambdap  <  0.0) 

•laa 


/••••  stable  condition  ••••/ 

(  aignac  «  1.0:  aignq  "  1.0;  I 
i  aignac  •  -1.0;  signq  -  -1.0;  ( 

(  signacp  ••  1.0:  signqp  -  ).0;  ) 

(  signacp  <■  -1.0;  signqp  >  -1.0;  I 

/•*••  unstable  condition  **•* 

(  Signac  >  -1.0;  signq  >  1.0;  ) 

I  Signac  •  1.0;  signq  -  -1.0;  ) 

(  signacp  >  -1.0;  signqp  •  1.0;  ) 

(  signacp  •  1.0;  signqp  •  -1.0;  I 


Ccfinal  •  (Eef  •<•0.0)  7  Ec(kount)  :  Ccf; 
Eepfinal  •  (Ecfp  ••  0.0)  ?  EcpIkount)  :  Ecfp; 


for (J«l; j<«kount; !♦♦) 

( 

Ech)rpa(j]  •  Ec(l]*Ecfinal/(Ec(l)«’(Ecfinal-Ecil])* 

•xp(acfactor*Ecfin-tl*xpI  j|) ) 

Ee)>ypadal( j)  -  Eehypal j]-Echypa(l); 

Eephypalj]  -  Ccp(l)*Ecpfinal/(Ecp(l)‘*'(Ecpfinal-Ecp|)  I )  * 

•xp(acpfactor*Ec|:[in.-<l*.';p|  j))  > 

Ecphypadall j]  -  Eephypal j]-Ecphypa(l]; 

Q  •  2.0*Eehypa( j)*Xstotal*(1.0-(lstotal/lttotal) ) 

-hsq*(Xstotal/Xttotal)*(1.0-(Xstotal/Xttr>t.tl) ) 

*  (latotal/Xttotal)  'SOdsp)  *50 (sigma)  ; 

Op  -  2.0*Eephypa|})*Xstotal*(1.0-(Xstotal/lttotnD ) 

-hsq*  (Xstotal/Xttotal)  •  (1 . 0-  ( Xstota  1/  It. t oi -i  I ) ) 
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ar« 

>  ((Ittecal/dttoCal-lateta))) 

*  ( dap*al«M)« (ai9nq‘sc|ri  (01 ) ) )  /  ('n' '  Hi-  -(> 

18  (axg  >  0.99999999)  arp-l.O; 

ataRO) 

•  57.29S779Sl31*acoa(aro) ; 

azfp 

*  ((Xttotal/dttotal-latetal)) 

*((-Ia*aioaa|*(Bi«n4p*s<itt  (Op)  I ) )  /  ('i-r  •  <<*•  •>' ' 
If  (ar«>  >  0.99999999)  acOp-l.O; 

atahpCjl 

*  S7.29S779S13i*acoa(axGP) : 

Ceaapljl 

*  (CelH-Cefinal)  *aap(aefactor*sp|  j) )  i-Rcr  i n.-i  i 

EcaapdaKjl 

-  EeaxplJl-Ceaapd); 

Cepaapt j) 

“  (Ecp|l)-Cepfinal)*aap(acpfactot*xp|  J))  in.-ii; 

EcpaapdaK]: 

-  Ccpai^ljl-Ccfwapll); 

0 

•-  2.0”Scaap(l)*Iatotal”d.0-d8total/lttotAl}) 

-haq* dacotal/Zttoeal) •d.O-datotal/lttot-i  1 )  I 
”  (Zatotal/lttetal)  *S0(Zap)  •SO(ai9M) 

0? 

-  2.0*Eepaap(J)*Zacecal*d.0-(latotal/utoini)  > 

-haq*  datetal/Ittotal)  *  d .  0-  (Istotal/ 1 1  t»>i  i )  i 
(latotal/Ittetal)  *$0(18)  *S0(siq«Ka)  ; 

az« 

•  ((Ittoeal/dttatal-latatal)) 

*((Z8p*8i9a«)4(8ionq*8qrt(Q))))/(»')ti  (h-t'})) 
it  (ac«  >  0.99999999)  arq-l.O; 

asaljl 

•  5Z,295779S13l*aeo8(ar9); 

•xgp 

-  ((Zttotal/(Ztcet8l-Z8total)) 

*((*l«*«i9M)i-(8t9nqp*sqct(0p))))/(!sqri.(li;i'|i ) 
if  <ar«p  >  0.99999999)  8r«p-1.0; 

atapi j) 

) 

”  57.29S779513l*aco8 (atgp) ; 

initiallta  graphics  •••••«•••••/ 

inlt  ( ) :  coJlor_aealt '  'cyanblua")  ; 

•calaCgraysealf :-) ;  wlndowOO ; 
b9Col(7);  axaaaO;  zalerCO); 

«ova(220,-:i0).-  pzlnz:  rHXHQOKX  OOAl  SPINNER  MASS  SPRING  SYSTEJl") 
xindow (22.  -S5. 979, 87;  ;  bgcol  (3) 

acala<0,10;;3. 10000.01;  zact(0.0.10000.10000l ; 

•  !•#  1 J  S 

.'3va(20,U0'»20Q);  data(); 
r5v«(-l«0,110«70|;  prlntf ("0"); 
ravadOfllO-OO);  printf(’'w3  -  w3  initial".'.- 

xaecorCO. IS.llO,  19);  rsva (-180, 19'»70);  princCCO")  ; 

nova  (10, 19-»200)  ;  pxiatf  (-xP")  ; 

vactor(0.1<.110.18);  .''3va(-180.  li^-YO)  ;  princf  CO”); 

r^ova  (10, 18«200) ;  pcintf  {”xl«)  ; 
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««ctwlt.  I>.ll0.t7);  pKiAtfl'O”!; 

■O««{10,17*2M| ;  prl*t((*'s  priaw  doi*); 

printfro*): 

■wv«(10.1(«200);  pciiiC(l*B  priM't; 

v««ter  (0.19,110,13):  «ov«(-lC0,  ISOO) :  prlMtf  (*0**) : 

■ov« (10, 19*200):  pclutCCt  doc"); 

vo«co'r(0,11, 110,14):  «ev« (-ICO,  14*70):  pciAte(*0*'): 

a»w«(10.14«200);  pclntC ("s") : 

««eter(0,13.110,13):  aevo  (-ICO,  13*70) :  printCCO"); 

Mvo  (10, 13*200):  prlncf  (")i«-ko  1  (r),  t>-r  i  tK)"); 

v«ctoc(0.12.110.12);  mow«  (-ICO,  12*70) ;  pciotCCO")  : 

M«« (10, 12*200) :  pclntC ("Ce-Ce  1  (x),  Ccttype-Cchypc  (  (it)"): 

v«etoc(0,ll,110.11):  mo«« (-ICO, 11*70) ;  pcintC(*0"): 

MOV* (10, 11*200) :  pxlatC ("Cep-Cep  1  (r),  Cep)iyp«-Ccphypc  i  O'l"). 

veetor (0,0, 110,0):  aeve (-ICO, 0*70):  pclntf("0"): 

Mvo (10,0*200);  pxlntf  ("tiMta,  eta  (b),  atap  (r),  )t  -  li  tnii 


vactex (0, 19*790, 11, 19*790)  : 
vactex (0, 19*900, 11. 19*900) : 
aactex (0, 19*290, 11. 19*290) : 
vaetex (0. 19-290, 11. 19-290) ; 
vaetox (0, 19-900, 11, 19-900) ; 
aaetox (0, 10*290. 11, 11*290) ; 
vaetox (0, 19-290, 11, 11-290) ; 
vactex (0, 19-300, 11, 19-900) : 

xseala  ••  10000. 0/xpOcoant); 


vaetox (19,19*750, 110, 19*750) 
vaetox  (19,19*500.  no,  19*500) 
vaetox (19, 19*250, 110, 19*250) 
vaetex(19,19-250,  no,  19-250) 
vaetex (19, 19-500, 110, 19-500) 
vaetox  (19. 19*290.  no,  19*290) 
vaetex (19, 19-250, 110, 19-290) 
vaetox (19, li-900, 110.19-900) 


eolex(l):  /•••««  piua  -  total  anaxoy,  xetex  )»ypa  •*•••/ 

Cox  ( )-l : JOtoont : J**) 

( 

voetex( (Int) (xaeala*i9( j] ) , (Int) (13*(0.01*atocaldal | j J ) ) , 

(Int) (xsealo*xp|)*l]) , (int) (13*(0.01*atotaldel()«l) > ) ) ; 
vaetox ( (Int)  (xscala*icp(  j) ) ,  (int)  (13*(0. l*atotaldal (}) ) ) . 

(int)  (xscala*i(p(  j*l]),  (int)  (13*(0.1*aCotaldel|  j*!))  > )  ; 

vaetox ( (int)  (xseala*xp(  jl ) ,  (int)  (12*(90.0*Cc>)ypadal  1  j) ) ) . 

(int)  (xscala*xp(  j*l|),  (int)  (12*(90.0*Cc)«ypedcl|j*>])>); 

vector ( (int)  (xscala-xptJJ).  (int)  (0*(100,0*ata)»|J))). 

(int)  (xacala«iip(]*l]).  (int)  (0*(100.0*atah(]*ll))); 
vaetox ( (int) (xacala'xpl]] ) , (int) (n*(1000.0*atah | j) ) ) . 

(int)  (xscala*xp(  j*n).  (int)  (n*(1000.0*ata)i|  j*l)) ) ) . 


color (5);  /•••••  pgjcpjy  «  platCoxm  exponential*****/ 

Cox  ( j-1; jOcount; j**) 

( 

vector ( (int) (xseale*xp( ]) ) , (int) (n*(50.0*Ccpaxpdel I jj ) ) . 

(int) (x8cala*xp( j*ll), (int) (n*(90.0*Ccpaxpdnlt  ) ) ) > ) ; 

vactor((int) (xxeala*xp|])), (int) (0*(100.0*atap( j))), 

(int)  (xacala*xp|]*l]),  (int)  (0*(100.0*etapM<l]) ) )  .- 
vector ((int) (xseale*xp(})). (int) (n*(1000.0*atap(}))). 

(int) (xxcala*xp(j*ll), (int) (ll*(1000.0*ctap| j»l ) ) | ) . 
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ceXextO;  /•••••  y«iiew  -  rotor  •■ponontlol*****/ 

tax  IJ-1; Jtkooot; !♦-» 

voetoxfUnt!  ':xsc«l«*xp(  jl).  tint)  (12«  ISO.O^Ecoxpdol  I  Jll). 

Cint!  (xscaio*xp(  Unt)  (12»<50.0*Ecexp«t«*l  I  i*  1 1 1  •  t . 

voctor ((inti ;xscolo*xp| j)) . (int ) (0« (100 .O*ota I j 1 ) 1 . 

•  (int‘  (xacalo*xp(l'»l|) .  (inti  (0*(l00.0*«ta|  j«-l  1 1 1  • 
voctox((int1 ixacalo*xp( j] ) . (inti (li«(l000.0*ctal  jl) >. 

(inf (xacala*xp| )«!)). (inti (ll*(1000.0*etai j* 1 ) M  » 


eoXox(0i;-  /•••••  black  -  piattora  hypotboiic  ••*«•/ 

tax  (J-1; J<kount; J**J 
I 

voctor  ( (inti  (xacaio*xp(  jl ) ,  (int )  (ll*(50.0*Ecphypcd<jl  I  jj ) ) . 

(in;;  (xscaXo*xpl  j«l]l  <  (inti  (il*  (50.0*Ecpliypcd«l  I  j*  1 1  n  )  ; 

voctor ( (int) <KscaXo*Kp{ j) I ,  (int) (0«(X00.0*otahp( j) ) ) , 

(int) (xacalo«xp(l«il) , (int) (0«(100.0*etahpl  j* 1 ) ) ) ) . 
voctor  ( (int)  (xaeaXo*xp(  j) ) ,  (int)  (li'*'(i000.0*otal«pl  jl ) )  > 

(int)  (xacalo*xp|  j^H).  (int)  (11«  (1000.0*eta)ip|  j»)  I ) ) ) . 


color (4);  /•••..  ^od  -  rotor  **•**/ 

tax  ( j-1; j<kount; j**) 

( 

voctor  ( (int)  (xaealo*xp(  jl),  (int)  (llO-t  (1432.39448783*  (yp(3)  I  jl  *vr !  -^11  > ) ) ) 
(int)  (xacalo*xp(  (int)  (110« (1432.39448783* (yp|3)  1  j*l) -ypl  3!  ( 1 1  >)  * 

voctor  ((int)  (xacalo*xp(  j)).  (int)  (19'*(1432.39448783*yp(21 1  j) ) ) . 

(int)  (xaealo*xp(  j«l)),  (int)  (19*-(i432.39448783*yp(2n  j • )  i ) ) ) 

vector  ((int)  (xacalo*xp(  jl),  (int)  (18'»(1432.39448783*yp|l)  t  j) )  > , 

(int)  {xacalo*xp(j4>ll),  (int)  (lB4(i432.39448783‘yp1 1 )  I  i O  ! ) ) ) 

voctor((int)  <x8calo*xp|  jl),  (int)  (174’(1000*yp{7)  i  jl) ), 

(int)  ixacalo*xp(  j*!)) ,  (int)  (n+dOOO^yplTl  ( j*  1)  1 1 1  : 

voctor(  (int)  <.'cacalo*xp(  j) ) ,  (int)  (16*(1000*yp|6)  I  jl) ), 

(int)  ;xacalo*xp(  j*'l)) ,  (int)  (IS*  (1000*ypj8i  I  j*l) ) ) ) : 

-^r{(int)  (xacalo*xp(  j)),  (int)  (15*  (1000*ypI5)  I  j))), 

(int)  (xacalo*xp(  j*!)) ,  (int)  (lS*'(1000*yp|S)  I  j*l  1 ) )  >  ; 

voctor ( (int) <xscalo*xp( j] ) , (int) (I4*(1000*ypl4 J  I  jl ) ) , 

(int;  ;xxcalo*xp(  j'*!)) ,  (int)  (14^(1000*yp|41  (  j' ) ) ) ) )  ; 

voctor  ( (int)  (xacalo*xp|  j) ) ,  (int)  (13'*(0.01*kedcll  j) ) ) , 

(int)  ;xacalo*xp(  j«l)) ,  (int)  (13'»  (0.01*kcdel )  j* )  1 ) ) )  ; 

vector  ( (int)  (xscalo*xp(  jl ) .  (int)  (13'*(0.1*kodel(  jl ) ) , 

(int)  ;xacalo*xp|  j-*l)),  (int)  (13'«'(0.1*kedell  j*!  | ) ) )  : 

vector! (int)  :xscalo*xp{  jl ) ,  (int)  (12-*(S0.0*Ecdol|  j)) ), 

(int)  ;xaealo*xp|]'fl)) ,  (int)  (12«(50.0*Ecdol(  j*  1 J ) ) )  : 

voctor  ((int)  (xscaXo*xp(  jl),  (int)  (ll-*(S0.0*Ecpdel(  j))), 

(int)  :x8ealo*xp( j*l|) , (int) (ll^iSO.O'Ecpdel I j<l)) ) > : 

vector  ( (int)  (xsealo*xp(  j] ) .  (int)  (0*-(1000.0*hdel  ( j ) ) ) . 

(int)  -'xaealo*xp(  j«l) ) ,  (int)  (0*'(1000.0*lidcl  I  j<  ))))): 
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wacteci  (inti  (Mcala'n^lll  I.  (Inc I  (0<  (iOO.O'Lbetal  J|l ) . 

(inti  (iiae«l«*npl)0  n .  (l«t)  (0»(XOO.O*th«t*U'  1 1 )»» 
v«ctec( (inti (K*eal«*npf  JJ I, (inti (11< (I000.0*thcra| jl 1 1 . 

(inti (xscnl«*npll«l|). (int) (ll*(i000.0*i  hetal j* 1 1 1 11 . 

/•  v«etoc( (int) (ase«l«*np(ll ) > (inti (0« (0.0001*ka( j|l I . 

lint)  (nseala'npt  (int)  l0*(0.000l*kel  J«l )))>. 

•/ 

/•  vactor ( lint) (xacnla'npl J) ). (int) (0*(0.0001*etotal I jl ) I . 

(int) (xaeala*npll«ll), (int) (0» (0.0001 'etotnl (}•) I ))) . 

*/ 

/•  vaetoc((int)  (Kacala*xp|  jl),  (int)  (0«(0.1*)t( ))) ) , 

(int) (xaeala*xp(j*l)). (int) (0*(0.1*h|)*l)>>); 

•/ 

) 

MindowOO : 
color (0) ; 

print  initial,  final  data  *“/ 


awvo(3SO.I7S);  pxintfl*  IHITIAL  COHOITIOMS  kODiii  -  tel* .  I.'mhh  I 


/• 

a)ova(22,910| :  pcintfCtiaw  •tS.Of  aaeonda'*.xptkount)): 

awvaiSOO.  910) ;  pcintf("apa  ••  to  %d/td'*.apt,nbad,  nok)  ; 

awva(22.940) : 

printfCHl  -%i.2f  ¥2  -94. 2f  113  -%4.2£",  yplDU) .  yp(2H  1 1 .  yi'l  M  i  1 1 )  • 
a»va (300,940);  printfl'wSp-  t9.7e  L  -  %3.1t".H3p,L) . 

aM>va(22,990):  printf("21  •%4.0f  12  •%4.0f  13  -tt.OI'*,  II.  12, )  i)  ; 

nova (22, 1010);  printf("M  -  94. Of  n  -97.Sf^M,m); 
nova (22, 1040);  printfCa  -94. If  k  -94. If  c  -94.1f'‘,a.k,c).- 

nova(22, 1070) ;  printfl"!  -94.2f  tdot-94.2f”,yp(4)  ll|  ,yp(Sni  1 1  •' 

nova (300, 910):  ptintf  ("Xlp-94.0f  12p-94.0f  13p»94 .Of*. Dp. I2r>  - 

nova  (300, 1010) ;  printf("Hp  -94. Of  ; 

nova (300, 1040);  printf  (*ap-94.1f  kp-94.1f  cp-94.1f*',ap,kp,et>> ; 

nova  (300, 1070);  printf(*xp  -94. 2f  tdocp -94.2f'',yp(0!  (n,yp(7)  (1)1  ; 

nova(22,1110);  printfCh  1-99. 4f  ka  i-97.1f’',))|  l),ka|l) )  ; 
nova (22, 1133);  printfCh  f-99.4f  ka  f-97.1f",h(kount),ke|kounc)) . 
nova (22, IISO) ;  printf ("hdalpar-9S .22", 100.0* (h (kount ] -h I 1 ) ) /h | 1 ) ) : 
nova (230, 1160) ;  printf  ("kadalpac-9S.2f",  100.0*  (ka(kount  l-ke(l  I )  /k»|  1 1 )  .■ 
lllOva  (300, 1160)  ; 

printf ("a  dalpac-9S.2f'',100.0*(atotal(kount]-atotal 1 1)) /"total  I ) 1 ) : 
nova  (730, 1160);  printf  ("Cep  dalpax-9S.2f",  100.0*  (Cep [ kount ) -Crp | )  | )  /Frpl  I ) ) 

nova (500, 1110);  printf ("Cc  i-99.4f  thaca  i-97.Sf".Cc(l),chata(] | ) . 

.isova  (300, 1133) ;  printf("Ec  f-99.4f  titaca  f-97.  Sf",Cc  I  kount],  theta  (^ouni  | )  .- 

•/ 


neva(22,900) ;  printf ("  platforn  rotor"); 

nova(22,920);  printfCXlp  •  96.2f",np); 

nova (240, 920) ;  printf (”Z1  -  96.2f",Zl); 

nova(22,940);  pxintf("Z3p  -  96.2C".I3p); 

nova (240, 940);  printf ("Z3  -  96.2f",I3); 

nova{22,960) ;  printf ("Hp  -  96.2f",Mp); 

nova  (240, 960) ;  printf(”H  -tO.Zf",!*); 

nova(22,990);  printf("np  -  96.2f",np); 

aova(240, 900) ;  printf("n  -  t6.2f",n); 

sova(22,1000);  printf ("ap  -  96.2f*,ap); 

nova (240. 1000) ;  ptintf ("a  -  96.2f",a); 

nova (22, 1020) ;  printf ("kp  -  96.2f",kp); 

nova (240, 1020);  prlncf("k  -  96.2f",k); 
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>v«(22. 1040) :  prlntfCep 
(240,1040):  pclnttCc 


-  16.2f",cp); 

-  \4.2f“,e); 


l■QV•(^2, 1060)  ; 
Mv«  (240, 1060). 


printfCL  -  «6.4C'‘.L); 
pclntecis/lt  -  \6.3I",  Istotal/Ittotwl) 


MOV* (22, 1100) ; 
MOV* (240, 1100) : 
Mov•(2^,11^0); 
Mov*'(240,1120): 
MOV* (22, 1140) 
Mov*(240, 1140): 
MOV* (22, 1160) ; 
MOV* (240, 1160); 


prlntt ("tp 
printt  Cs 
printt ("dtp 
pcintf  Cdt 
print ( ("wSp 
princtCwS 
print! ("wl 
print!  (''w2 


%6.2I-,ypl6)|l)); 
t6.2I‘,yp(4| (l|); 
%6.2!-,ypn)Cl)); 
16.2!-,yplS)U)); 
«6.2f*‘.yp|3)  (l)«si9iKt); 
%6.2f‘,yp(3)il)); 
\6.2('‘,yp(l)|l)): 
t6.2f-.ypl2)|l)); 


MOV* (400, 900); 
MOV* (130, 900) : 
MOV* (480, 920) ; 
MOV* (130, 920) : 
MOV* (480, 940); 
MOV* (130, 940) ; 
MOV* (480, 960) ; 
MOV* (130, 960) : 
MOV* (480, 980) ; 
MOV* (130. 980) ; 
1(480,1000): 
»v*(130, 1000) : 
'print!  (' 


print!  (''%3.0!  s*conds'*,  nplkount ) ) . 

print!  Ch  i  .  -  -  — . 

print! ('id  /  %d  / 
print!  Ch  ! 
print!  C*p«  - 
print!  Chd*lp  - 
print  !C)c*  i 
print! ("*  i 
ptint!(’'k*  !  - 

print! '‘'*  !  - 

print  !  ("(cndnlp  - 


%9.4£-,h(l)); 

%d",  nb«d.  no)c.  kount) : 

«9.4!*',h(kount)): 

%1*'‘,*pt); 

%5.3!*.100.0*(h(kount  l-)«|1 1  )/hi  '  ! . 

kf.ii'.k*!!)); 

tO-ii^.^totaldl): 

«9.4!‘',k*Ikount]); 
\9.4!'',*totaltkoant))  ; 
»5.3!",100.0»()t*Ikount)-ke()  1 )/!'')  I  i). 


MOV* (480, 
MOV* (130, 
MOV* (480, 
MOV* (130, 
MOV* (480, 
MOV* (130, 
MOV* (480, 
MOV* (730, 
MOV* (480, 
MOV* (730, 


1030) ; 
1030) ; 
1050); 
1050) : 
1070) ; 
1070) ; 
1090) : 
1090) ; 
1110); 
1110); 


*«d«lp  -  OS.II”,  100.0*  (*tot«l(kount)-«totAl  1 1 1 ) 1 1  I ) 

print! CEcpl  -  %9.3f‘‘,Ecp(l)); 
printfCEc  i  -  %9.3£".Ec(l)) : 
print!CEcp!  -  %9.3!*',Ecp(kount) ) : 
print! (”Ee  !  •  %9.3I",Ec|)(Ount)) ; 

print !(-l«nibd«p  •  %5.3!*’,l«ii8)d«p); 
print!ClaMbd«  -  «5.3!‘‘, laMbd*) ; 
print! Cslfnqp  -  %5.2!'',si9nqp); 
print!("si9nq  -  95.2f",si9nq); 

print! CdSep/lambdap*  %6.3!",  (Ecp()toont)-Ecpr)))/l.''n!'.iip)  ; 
print! ("dEc/ lambda  -  96.3!",  (Ec(ko«nt)-EcIl  J ) /ln»r*''l.i) . 


MOV* (480, 1130);  print! (’*Ee!p  -  %8.6!",Ec!p); 

MOV* (730, 1130);  print! CEc!  -  i8.6!‘',Ec!); 


Mov* (480, 1155);  print!Cth*t*  i  -%6.2£'’,th«t*|l))  ; 
mov*(620,1155);  print!C*tap  i  -96.2!’',*tap(l])  ; 

MOV* (760, 1155);  print!C«ta  i  -t6.2!”,*ta|l]); 

MOV* (480, 1175):  pcint!Cth«ta  !  •96.2t'',th*ta(kount])  : 
MOV* (620, 1175);  print!C*tap  !  •t6.2!”,*tap(kount)); 
MOV* (760, 1115);  print!C*ta  !  ••%6.2!",*ta|kount)); 


/•••  upp*r  rlqht  hand  cornar  •••/ 

MOV* (600, >33);  printiCEc  initial/!inal  -  kq  /  %9",EcllJ,Ec(kount; ) ; 
MOV* (600, -13);  printiCEcp  initial/!inal  -  iq  /  9q”,Ecp|l),Ecp():cun;.  i )  ; 
mov*(600,27)  :  print!CEc  !actor  /  Ec!  -  %q  /  %q'',«c£actor,Ecf  >  ; 
Mov*(600, 47) ;  print!("Ecp  factor  /  Ec!p  -  %q  /  »9",ecp!actor.Ecfp) ; 

MOV* (980. -33) ;  print! Ctd'.nbad) : 
mov*(980,-13) ;  print !( "kd", nok ) ; 
mov*(980,1);  print!  ("kd’’,  kount  ); 

MOV* (980, 27);  print! (“nbad**); 
mov*(980, 47);  print! ("nok") ; 

MOV* (980, 67) ;  print! Ckount") ; 

MOV* (980, 97);  print! (-96. 4!", (100.0* (kelkount) -kis(l] ) ) /kc| ) | )  ; 
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Mw«(ftO. n't!;  print* ini)/-  '  iMHi 
WvclMO.  X3'>» ;  print*  (XOO.O*  Ihlkounl  l-h|  1 11 )  /hi  1 1 )  : 

mw«(9I0.1S7);  print*  (-kndvl'*); 

■ov«(M0.177);  print* ('•dnlp'); 

MOv«(9e0.X97);  print* rhdnXp-); 

n»v«  (910, 227) ;  print* ("kf". that* ( 1 1 ) : 

■ov«(9l0,247) ;  print*('*k9*‘,th«tn(kottnt)).* 

•ova (990, 297);  print* ("th  X-); 
aiova(990,297):  print*  Cth  *-); 

a»va(9IO, 3X7) ;  print* (-«9''.ata(X] ) ; 
niova(980.  337) ;  print*  ('‘%9'',ataikount| )  ; 

nova (980,357);  print* ("ata  1"); 
nova  (980, 377);  print*  ("ata  f); 

nova (980, 907) ;  print* (''%9",atap(l] )  ; 
nova (980.  427) ;  print*  (**%9*',atap(kottntl ) : 

nova (980. 447) ;  print*  ( "at apX ")  ; 
nova (980,  467) ;  print* ("atap*") ; 

iwva(980,497) ;  prlnt*("%9",Ce(koantl*Ec(X)) ; 
nova (980, 5X7);  print* ("%9”,Xaiibda); 
nova (980, 537) ;  print* (”%9", sl9nq) ; 

nova (980, 557);  prlnt*("Ecdal"); 
nova  (980, 577);  print*  ("XaaMa")  ; 
nova  (980,  597) ;  print*  ("sipnq") ; 

nova  (980,627);  pr Int* ("%9", Ecp ( keunt ) -Eep | X 1 ) ; 
nova (980, 647) ;  print* ("kq", Xanbdap) ; 
nova  (980, 667) ;  print* ("kq", aiqnqp) ; 

nova  (980, 687);  print*  ("EepdaX") ; 
nova  (980, 707);  print*  ("XanMap") ; 
nova (980, 727);  print* ("aiqnqp") ; 

nova (980, 757) ;  print* ("kq",  IstotaX/lttotaX) ; 
nova  (980, 777) ;  print* ("la/It") ; 

fraa^matrix  (yp,  X,K,  :,MAXAIU«AY) ;  *raa_vactor  (Ec,  1,ma.':appa') 

fraa^vactox  (xp,  X,IIAXAIUtAY) ;  *raa~vactor  (Ecdcl,  1,KAX At- >■•*•.?>  ; 

iraa^vactor  (yatart,  1,N) ;  Iraa^vactor  (Ecp,  1,HAXARPAV)  ; 

fcaa^vactoc  (lea,  l,HAXAfU(AY) ;  Ccae^vactor  (Ecpdel,  1,kaxaff.ay)  ; 

*raa~vaetoz  (kadal,  X,MAXARRAY) ;  *cee~vactor  (Echypa,  I, MAXAFRAY )  ; 

*raa~vactQz  (atotaX,  1,MAXARRAY) ;  iraa^vactor  (Ecliypedcl .  i.maxaf.i'ay)  : 

2zaa~vactoz  (atotaXdal,  X.HAXARRAY) :  iraa^vactor  (Eephypo,  1,haxa".i>ay»  ; 

*zaa~vactor  (h,  X.HAXAFRAY) ;  *caa~vactor  (Ecphypcdal,  x,r.AXARRAV) 

*zaa~vactoz  (thata,  1,maxarray)  ;  fraa^vactoz  (eta,  1,haxapray)  ; 

*zaa~vaetor  (hdaX,  1,  KAXAIWAY) ;  fraa^vactor  (atap,  X ,  haxarkay ) 

izaa'vactoz  (Ecaxp, X,XAXA1UIAY) ;  *cac~vactor  (Ecaxpdel,  x,:ia>:afi'AY)  ; 

*raa~vectot  (Ecpaxp,  l.MAXARRAY) ;  traa~vactoc (EcpexpdQl ,  l.tiAXAKi’AY) 

*zaa*'vaetoc  (atah,  l,KAXAlUUiY) ;  *raa” vector  (ctahp,  l,HAXAi'.i!AY) ; 
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